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ABSTRACT 


A  magnetic  fuel  cell  uses  high-pressure  gas  to  cause  a  time-dependant  change  in 
magnetic  field  to  generate  a  voltage  based  upon  Faraday’s  Induction  Law.  In  this  thesis, 
the  use  of  a  compressed  magnetic  gas  passing  through  a  magnetic  field  is  described  for 
generation  of  electrical  power. 

The  first  stage  of  this  thesis  explores  the  fundamental  physical  nature  of 
diamagnetic  and  paramagnetic  gases.  Once  this  stage  is  completed,  it  will  become  clear 
that  there  are  several  methods  that  can  be  used  to  overcome  the  weak  nature  of  these 
gasses  and  achieve  useful  work,  due  to  pressure-sensitive  changes  in  magnetic 
susceptibility.  The  second  stage  of  the  thesis  develops  the  theory  used  to  model  the 
operation  and  physical  characteristics  of  the  magnetic  fuel  cell.  The  third  stage  shows  the 
development  of  the  specialized  magnetic  fields  required  and  equipment  configuration 
necessary  to  employ  the  above  methods  and  predicts  the  experimental  measurements  of 
the  resulting  voltages  and  currents.  The  fourth  stage  discusses  experimental  findings  and 
presents  an  optimized  configuration.  The  final  stage  of  the  thesis  concludes  with  lessons 
learned  and  interesting  ramifications  of  the  fuel  cell. 
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EXECUTIVE  SUMMARY 


A  magnetic  fuel  cell  uses  high-pressure  gas  to  cause  a  time-dependant  change  in 
magnetic  field  to  generate  a  voltage  based  upon  Earaday’s  Induction  Eaw.  In  this  thesis, 
the  use  of  a  compressed  magnetic  gas  passing  through  a  magnetic  field  is  described  for 
generation  of  electrical  power. 

The  first  stage  of  this  thesis  explores  the  fundamental  physical  nature  of 
diamagnetic  and  paramagnetic  gases.  The  ability  of  these  gases  to  change  their 
susceptibility  as  pressure  changes  was  instrumental  in  producing  voltage  in  the  fuel  cell. 
The  basic  equation  explaining  how  susceptibility  affects  the  magnetic  field  of  a  magnetic 
circuit  is  show  in  Equation  (1.1), 

b  =  mS^  +  x)h.  (1.1) 

The  magnetic  flux  density,  B  ,  equals  a  constant  ( ju^  which  is  the  permeability  of  free 
space)  times  the  quantity  one  plus  the  susceptibility,  x  >  the  magnetic  field  strength, 
H .  This  equation  explains  how  a  change  in  susceptibility  can  affect  the  magnetic  flux 
density  of  the  magnetic  circuit.  The  final  component  of  the  magnetic  fuel  cell  concept  is 
gas  dynamics.  Using  Bernoulli’s  principle,  the  velocity  of  a  high-pressure  gas  can  be 
determined.  Velocity  of  the  high  pressure  gas  and  Equation  (1.1)  can  be  combined 
together  and  substituted  into  Earaday’s  Induction  law  to  achieve  a  linear  relationship 
between  the  voltage  output  of  the  fuel  cell  and  the  gas  pressure. 

Once  the  basic  derivation  was  completed,  it  is  clear  that  there  are  several  methods 
that  can  be  used  to  overcome  the  weak  nature  of  diamagnetic  and  paramagnetic  gases  and 
achieve  useful  work.  The  second  chapter  of  the  thesis  developed  the  theory  used  to 
model  the  operation  and  physical  characteristics  of  the  magnetic  fuel  cell.  The  predicted 
voltage  output  was  3.6  millivolts  with  a  period  of  3.7  milliseconds.  A  schematic  of  the 
test  setup  is  shown  in  Eigure  1.  The  largest  challenge  during  this  testing  phase  was 
reducing  the  noise  floor  and  accounting  for  mechanical  vibration  produced  by  the 
discharge  of  the  high-pressure  gas. 
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The  schematic  (Figure  1)  shows  the  gas  flow  path  for  the  high-pressure  air.  The  gas 
flows  past  the  pressure  regulator  through  valves  one  and  two  into  the  magnetic  fuel  cell. 
The  gas  inside  the  fuel  cell  modulates  magnetic  flux  inducing  a  voltage  in  the  wire  coil. 
The  oscilloscope  records  the  voltage  signal.  The  microphone  triggers  the  oscilloscope  to 
start  recording  the  signal.  The  isolated  voltage  signal  of  6.3  millivolts  produced  by 
discharging  1800  psi  air  out  of  the  magnetic  fuel  cell  as  shown  in  Figure  2.  The 
sinusoidal  pulse  had  a  period  of  0.3  milliseconds. 

File  Control  Setup  Measure  Analyze  Utilities  Help  5:07  PM 


Figure  2  TIP  Oscilloscope  Screen  Capture  of  Magnetic  Fuel  Cell  Output. 


XVI 


The  vertical  scale  was  20  millivolts  per  division  and  the  horizontal  scale  was  100 
microseconds  per  division.  The  difference  between  theory  and  measured  values  was  3.6 
millivolts,  an  error  factor  of  0.57.  The  predicted  and  measured  voltage  showed  good 
agreement.  The  period  of  the  actual  signal  was  off  by  a  factor  of  a  hundred  but  agreed 
more  closely  with  the  manufactures  calculations  of  0.5  milliseconds. 

Chapter  V  summarized  the  first  three  chapters  and  discussed  a  recommendation 
for  the  next  experiment.  Chapter  V  concluded  with  ramifications  of  the  magnetic  fuel 
cell.  The  magnetic  fuel  cell  could  be  optimized  to  throttle  and  stop  the  high-pressure  gas 
flow  without  physical  contact.  The  equations  for  the  magnetic  fuel  cell  also  could  be 
applied  to  understanding  weather;  explaining  weather  effects  such  as  lightening,  storms, 
how  wind  is  generated  and  how  the  ionosphere  becomes  electrically  charged. 
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I.  INTRODUCTION 


Modern  portable  energy  generation  methods  typieally  rely  on  ehemieal  energy  to 
produee  eleetrieal  energy.  Fuels,  sueh  as  hydrogen,  oxygen,  diesel,  and  gasoline,  are 
used  to  generate  power  for  vehieles  and  other  types  of  portable  equipment.  These 
methods  bring  several  disadvantages;  they  are  highly  eombustionable,  explosive,  toxie, 
and  eostly  to  produee.  All  generate  environmentally  unsafe  waste  produets  (e.g.,  CO, 
CO2,  NO2  ete.). 

Portable  energy  generation  systems  based  upon  ehemieal  energy  are  very 
ineffieient.  After  the  ehemieal-based  fuels  release  their  energy,  the  energy  must  go 
through  several  eonversions  before  being  turned  into  eleetrieal  energy.  For  example,  a 
turbine  generator,  found  on  many  naval  ships,  has  very  low  effieieney.  The  turbine 
generator  proeess  eonverts  gasoline  into  meehanieal  energy  by  ignition  (with  some 
energy  lost  due  to  heat).  The  resulting  meehanieal  energy  is  then  ehanneled  through 
turbine  wheels  where  some  of  the  energy  is  lost  to  frietion  and  windage  losses.  The 
turbine  axle  rotates  an  eleetrieal  generator  that  has  frietion  losses,  hysteresis  losses,  and 
eounter-induetion  losses.  All  these  losses  oeeur  before  useful  eleetrieal  energy  is 
developed.  The  total  effieieney  of  the  system  ean  be  found  by  multiplying  individual  loss 
effieieney  terms  together.  The  resulting  total  effieieney  term  for  the  turbine  generator  is 
very  low. 

The  optimal  heat  engine,  of  whieh  most  ehemieal  energy  systems  are  a 
subeategory,  also  has  a  limit.  This  limit  is  based  upon  a  eoneept  ealled  Carnot  effieieney 
[1].  The  Carnot  effieieney  is  found  by  dividing  the  temperature  differenee  between  a 
heat  souree  and  a  heat  sink  by  the  temperature  of  the  heat  souree.  Using  a  praetieal  heat 
engine’s  temperature  parameters,  a  heat  engine  eould  aehieve  a  Carnot  effieieney  of 
approximately  45%  [1]. 

The  effieieney  limitations  of  the  heat  engine  are  one  of  the  main  reasons  eleetrie 
fuel  eells  have  been  mandated  by  President  Bush’s  administration  as  the  teehnology  of 
the  future  [2].  A  fuel  eell  has  a  one-step  eonversion  proeess  that  is  not  based  upon  a  heat 
eyele  and  therefore  does  not  have  the  Carnot  effieieney  limitation.  Unfortunately,  most 
eleetrie  fuel  eells  use  hydrogen,  whieh  is  explosive.  This  has  affeeted  the  growth  of  this 
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technology  in  the  commercial  sector.  The  main  research  effort  in  industry  has  been  the 
search  for  a  safe,  non-explosive,  low-cost  hydrogen  storage  tank  to  circumvent  this  issue. 
(Note:  Appendix  A  discusses  hydrogen-based  electric  fuel  cells  in  detail.) 

Instead  of  using  chemical  energy  in  a  fuel  cell,  it  is  possible  to  use  a  magnetic  gas. 
Direct  conversion  from  mechanical  energy  into  electrical  energy  is  achieved.  This  fuel 
cell  is  called  a  magnetic  fuel  cell.  Like  the  electric  fuel  cell,  it  is  not  a  heat  engine. 
Therefore,  the  theoretical  efficiency  could  be  much  higher. 

Magnetic  devices  and  applications  have  experienced  a  tremendous  growth 
commercially  in  the  last  forty  years  mainly  due  to  the  ability  of  industry  to  manufacture 
stronger  and  more  powerful  permanent  magnets.  In  the  early  nineteenth  century, 
permanent  magnets  were  hard  to  obtain.  Most  scientists  used  direct-current  transformers 
to  achieve  strong  magnetic  fields.  Permanent  magnets  were  only  several  hundred  Gauss 
strong  and  weakly  coercive.  (Coercivity  is  defined  as  a  property  of  magnetism,  which 
relates  to  how  resistant  the  magnet  is  to  demagnetization  by  an  external  magnetic  field.) 
Once  the  discovery  of  the  Neodymium-Iron-Boron  (abbreviated:  Nd-Fe-B)  and 
Samarium  Cobalt  (abbreviated:  Sm-Co)  rare  earth  magnets  were  made,  magnetic  flux 
densities  rose  within  a  few  years  to  10  to  20,000  Gauss  (I  to  2  Tesla).  The  discovery  of 
these  rare  earth  magnets  also  allowed  the  size  of  the  magnets  to  decrease  without  a 
decrease  in  strength  and  coercivity. 

To  achieve  voltage  induction  without  moving  parts  requires  the  development  of  a 
magnetic  gas  as  the  working  medium.  Magnetohydrodynamics  (MHD)  offers  one 
solution.  MHD  uses  ionized  gasses  and  liquids  as  the  working  medium  but  this  requires  a 
method  to  ionize  the  gasses/liquids  [3].  What  about  using  a  magnetic  gas?  The 
traditional  method  of  manufacturing  a  permanent  magnetic  material  involves  the 
establishment  of  magnetic  domains  within  the  material.  These  magnetic  domains  in  a  gas 
would  change  a  gas  into  either  a  liquid  or  solid.  Magnetic  liquids  have  been  made  using 
finely  powered  magnetic  materials  placed  in  an  oil  slurry.  These  liquids  are  employed 
effectively  today  as  new  generation  shock  absorbers  on  cars  and  as  hydraulic  lift 
replacements  for  carrier  launch  pads.  The  problem  with  a  magnetic  liquid  is  the 
difficulty  in  obtaining  it  and  it  is  usually  limited  in  quantity. 
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The  only  types  of  gas  that  exhibit  magnetie  properties  are  diamagnetic  (e.g., 
nitrogen)  and  paramagnetic  gases  (e.g.,  oxygen).  Diamagnetic  gases  are  gases  that  are 
repelled  from  magnetic  fields.  The  electrons  in  a  diamagnetic  material  rearrange  their 
orbits  slightly,  creating  small  persistent  currents  that  oppose  the  external  magnetic  field. 
Paramagnetic  gases  are  gases  that  are  attracted  towards  magnetic  fields.  Both  magnetic 
gases  are  weakly  magnetic  as  categorized  by  their  magnetic  susceptibility  (measure  of  a 
substance’s  ability  to  conduct  a  magnetic  field).  Paramagnetic  gases  are  much  stronger 
than  diamagnetic  gases. 

A  schematic  of  several  energy  transformation  concepts  is  shown  in  Figure  1.  This 
schematic  shows  that  useful  work  from  a  heat  engine  begins  as  chemical  energy  and  is 
converted  using  chemical  combustion.  The  combustion  energy  then  applies  mechanical 
force  to  move  a  prime  mover.  The  prime  mover’s  energy  is  then  electromagnetic  ally 
converted  to  useful  work. 


♦  Heat  Engines  -  Carnot  Efficiency 


♦  Hydrogen  Fuel  Cell 
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Work 
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♦  Magnetic  Fuel  Cell 


Figure  1  Energy  Transformation  Concepts. 


The  hydrogen  fuel  cell  produces  useful  work  by  converting  stored  hydrogen 
catalytically  into  electrons  and  protons.  The  electrons  flow  through  a  circuit  around  the 
fuel  cell  and  the  protons  diffuse  through  a  proton  exchange  membrane.  On  the  other  side 
of  the  proton  exchange  membrane,  the  protons  combine  with  air  to  produce  water  and 
heat.  The  electrons  in  the  circuit  produce  a  DC  current,  which  provides  work. 
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A  magnetic  fuel  cell  uses  the  varying  susceptibility  of  high-pressure  gases  as  its 
energy  source.  When  exposed  to  a  magnetic  field,  these  gasses  modulate  the  magnetic 
flux  density  and  induce  an  AC  voltage  in  a  winding.  The  AC  voltage  is  used  to  produce 
work. 

Finally,  a  battery  is  used  to  store  chemical  energy.  The  chemical  energy  is 
released  upon  closure  of  the  electrical  circuit  via  electrolytic  charge  transfer.  A  DC 
current  is  produced  which  is  used  for  work. 

With  an  understanding  of  the  energy  transformation  process.  Chapter  II  addresses 
how  the  varying  susceptibility  of  high-pressure  gases  produces  work.  Chapter  II  begins 
with  the  theoretical  explanation  of  magnetism  and  gas  dynamics  and  finishes  with  a 
presentation  of  the  magnetic  fuel  cell  concept.  Chapter  III  starts  with  a  full  theoretical 
development,  including  the  magnetic  circuit.  A  discussion  is  included  in  optimizing  the 
number  of  windings  of  the  magnetic  fuel  cell.  Chapter  III  ends  by  discussing 
implementation  considerations  of  the  experiment  and  experimental  predictions.  Chapter 
IV  addresses  the  experimental  setup  and  includes  an  analysis  of  the  data.  Chapter  IV 
concludes  with  a  discussion  of  the  lessons  learned  from  the  experiment.  Chapter  V 
provides  a  conclusion  with  several  recommendations  regarding  future  work  and  a  few 
ramifications  of  the  magnetic  fuel  cell  concept. 
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II.  BACKGROUND  THEORY 


This  Chapter  develops  the  underlying  physical  principles  of  how  a  magnetic  fuel 
cell  works.  The  first  section  discusses  magnetism.  Specifically,  it  describes  magnetic 
materials  and  their  classification,  the  concept  of  coercivity  and  susceptibility,  the 
hysteresis  loop,  and  magnetic  circuits.  The  next  section  discusses  gas  flow  dynamics, 
Bernoulli’s  principle,  the  continuity  equation  and  how  to  estimate  gas  velocities.  Section 
C  describes  Faraday’s  Induction  law,  showing  how  to  achieve  a  voltage  using  a  change  in 
the  magnetic  field.  The  final  section  introduces  the  basic  theoretical  idea  of  the  magnetic 
fuel  cell.  It  will  show  how  magnetism,  gas  dynamics  and  Faraday  Induction  concepts  are 
combined  together  to  achieve  the  overall  concept. 


A,  MAGNETISM 

One  fundamental  force  of  nature  is  magnetic  force,  a  property  of  magnetism. 
Magnetism  allows  the  application  of  force  at  a  distance  and  the  ability  to  connect  two 
metallic  objects  without  using  physical  fasteners.  To  understand  what  makes  a  material 
magnetic  requires  the  development  of  two  magnetic  concepts,  atomic  magnetism  and 
molecular  magnetism. 

Atomic  magnetism  focuses  on  the  structure  of  the  atom.  An  atom  is  made  up  of 
neutrons  and  protons  in  the  nucleus  and  electrons  orbiting  the  nucleus.  The  electrons 
have  a  property  called  spin.  This  property  has  been  identified  as  the  essential  ingredient 
to  producing  magnetic  fields.  Spin  can  be  broken  down  into  two  states,  spin  up  and  spin 
down.  An  elementary  rule  of  electron  spin  is  that  electrons  like  to  pair  up  with  one  spin- 
down  electron  and  one  spin-up  electron.  When  an  electron  does  not  pair  up,  the  element 
will  be  considered  magnetic.  Most  elements  are  not  magnetic.  Typically  unpaired 
electrons  are  in  the  atom’s  outer  orbit.  These  electrons  pair  up  with  electrons  from  other 
atoms  forming  non-magnetic  molecules.  The  exception  to  this  is  the  magnetic  elements, 
two  of  which  are  iron  and  nickel.  Iron,  for  example,  is  magnetic  because  there  are 
several  unpaired  electrons  in  its  3D  atomic  shell.  The  outer  shell  of  Iron  is  4S.  The  outer 
shell  acts  as  a  shield  for  the  D  shell  electrons  preventing  pairing  of  electrons  from  other 
atoms.  However,  not  all  iron  materials  are  magnetic  due  to  molecular  magnetism. 
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When  atoms  combine  into  molecules,  or  form  into  crystals,  only  certain  magnetic 
orientations  are  allowed.  Two  adjacent  atoms  with  unpaired  D  shell  electrons  will  have 
magnetic  fields  pointing  in  parallel  but  opposing  directions.  This  effectively  cancels  out 
the  net  magnetic  field.  Thermal  noise  constantly  overpowers  the  weak  magnetic  fields  of 
individual  atoms.  The  material  then  develops  regions  or  domains  within  it  that  are 
magnetically  polarized  in  different  directions  with  no  overall  magnetic  field.  Only 
certain  materials  exposed  to  an  external  magnetic  field  will  retain  a  residual  magnetic 
field  after  the  external  magnetic  field  is  removed. 

1.  Magnetic  Materials 
(a)  Classification 

The  following  categories  show  how  a  magnetic  material  is  classified. 

•  Ferromagnetic  materials  are  materials  in  which  magnetic  domains  align  to 
form  a  residual  magnetic  field.  This  magnetic  field  will  exist  without  an 
external  field  applied.  One  example  of  this  is  iron.  Iron  has  four  unpaired 
electrons  in  its  atomic  D  shell,  which  are  protected  by  an  outer  shell.  These 
electrons  have  the  same  spin  state  (up  or  down),  producing  a  strong  magnetic 
field. 

•  Anti-ferromagnetic  materials  are  materials  in  which  internal  molecular 
magnetic  fields  align  in  parallel  but  opposite  directions.  The  result  is  a  zero 
net  magnetic  field.  It  is  very  difficult  to  change  this  magnetic  field  using  an 
external  magnetic  source.  This  molecular  phenomenon  can  be  seen,  for 
example,  in  copper  atoms  in  a  metallic  lattice. 

•  Diamagnetic  materials  are  materials  that  have  no  internal  magnetic  field. 
When  an  external  magnetic  field  is  applied,  internal  magnetic  fields  are 
generated  to  repel  the  external  field.  The  atoms  within  these  materials  have 
no  unpaired  electron  spins.  Bismuth,  having  the  largest  number  of  paired 
electrons  in  the  periodic  table,  is  the  strongest  diamagnetic  material.  Striated 
carbon  has  been  reported  as  stronger  but  only  in  one  crystal  plane  [4]. 

•  Paramagnetic  materials  are  materials  where  a  majority  of  atoms  have 

magnetic  fields  aligned  in  different  directions  due  to  thermal  noise.  This 

6 


produces  a  zero  net  magnetic  field.  When  an  external  field  is  applied,  these 
internal  magnetic  fields  slightly  align.  Oxygen  is  one  element  that  is 
paramagnetie. 

•  Ferri-magnetic  materials  are  materials  where  an  external  magnetie  field  will 
cause  massive  alignment  of  internal  magnetic  fields  until  the  external 
magnetic  field  is  removed.  Elements  such  as  manganese,  aluminum  and 
cobalt  are  ferri-magnetic. 

(b)  Soft  or  Hard  Magnetic  Materials 

Magnetic  materials  are  elassified  as  soft  or  hard  depending  on  their 
eoercivity  ( ),  and  permeability  {ju).  Coereivity  was  defined  in  Chapter  I.  A  material 

with  a  high  coereivity  will  require  a  large  external  magnetie  field  to  reverse  its  magnetic 
polarization.  A  material  with  low  eoereivity  will  require  a  small  external  magnetie  field 
to  reverse  its  polarization.  Permeability  is  a  constant  of  proportionality  between 
magnetie  flux  density  and  magnetic  field  strength.  It  shows  how  effectively  a  material 
channels  a  magnetie  field.  Hard  materials  have  a  low  permeability  and  a  high  coereivity 
( >10000  A/m)  [5].  Soft  materials  have  a  high  permeability  and  low  coereivity 

(//^<1000  A/m)  [5].  For  the  magnetie  fuel  eell,  the  materials  of  interest  are  soft 

magnetie  materials.  The  magnetic  fuel  cell  is  designed  to  operate  in  the  kHz  range  with  a 
constantly  changing  magnetic  field.  Therefore,  materials  with  high  permeability,  but 
very  low  eoereivity,  are  desirable. 

2.  Susceptibility 

A  measure  of  a  material’s  ability  to  repel  an  external  magnetie  field  is  ealled 
susceptibility  ( j).  Susceptibility  is  unit-less  and  seen  in  Equations  (1.1)  and  (1.2);  as 

either  relating  to  relative  permeability  (//^)  or  a  constant  of  proportionality  between 
magnetie  flux  density  {B)  and  magnetic  field  strength  ( // ), 

A.=l  +  J,  (1-1) 
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B  =  MoMrH , 


(1.2) 


where  ju^  is  the  relative  permeability  of  free  spaee 


kg-m 
A' -s' 


=  4;r  X 10 


Experimental  susceptibility  (/r)  has  units  of  cm'/mol .  The  equation  for  experimental 
susceptibility  is  (MKS), 


=  X 


M 

Anp 


(1.3) 


where: 


•  M  is  the  molar  mass  (g/mol)  and 

•  /?  is  the  material  density. 

The  basic  equation  that  relates  magnetic  flux  density  ( 5  )  to  magnetic  field  strength  {H) 
and  magnetization  ( M  )  is: 

B  =  p{H  +  M).  (1.4) 

Susceptibility  is  related  to  the  basic  magnetic  equation  by  observing, 

B  =  (1-5) 

where  p  is  magnetic  permeability. 

The  form  used  herein  will  be: 

B^pS^  +  X)H.  (1.6) 


(a)  Diamagnetism 

In  diamagnetic  materials,  susceptibility  values  are  negative  and  very  small 
ix  ~  -10^^  [1]).  The  strongest  diamagnetic  material  is  striated  carbon,  x  ~  -0.00142  [4]. 
There  are  several  diamagnetic  gases,  a  few  examples  include:  fluorine,  steam,  and  the 
noble  gases.  Nitrogen  has  a  susceptibility  of  j  =  -5.4  x  10  [6]  at  room  temperature  and 
pressure.  Plasma  is  considered  weakly  diamagnetic  with  a  susceptibility  close  to  that  of 
nitrogen.  Since  diamagnetism  occurs  at  the  atomic  level,  quantum  mechanics  is  used  to 
explain  the  theory  of  diamagnetism.  Materials  become  more  diamagnetic  as  size  shrinks 
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towards  the  atomic  state.  Susceptibility  can  be  calculated  using  the  following  equation 

m, 


6m, 


where; 

•  is  atoms  per  unit  volume, 

•  ^  is  the  electron  charge  =  1.6  x  10“''* coulombs) , 

•  <  >  is  the  expectation  value  of  the  orbital  area  , 

•  Z  is  atomic  number  of  the  element  and 

•  m,  is  electron  mass  (w,  =  9.1  x  10  ^*kg) . 


(1.7) 


Susceptibility  is  therefore  proportional  to  atoms  per  unit  volume.  The  susceptibility 
equation  predicts  that,  by  keeping  volume/temperature  constant,  a  steadily  increasing 
pressure  will  cause  a  steadily  increasing  susceptibility.  Nitrogen  is  a  gas  that  makes  up 
approximately  80  percent  of  the  Earth’s  atmosphere.  Figure  2  shows  pressure  versus 
susceptibility  for  diamagnetic  materials:  nitrogen  (N2),  bromine  (Bri),  bismuth  (Bi);  and 
the  paramagnetic  gas,  oxygen  (O2).  As  pressure  and  susceptibility  increase,  bromine’s 
susceptibility  becomes  larger  than  bismuth.  Nitrogen,  on  the  other  hand,  must  be 
pressurized  to  several  million  psi  before  becoming  stronger  than  bismuth.  In  contrast  to 
these  diamagnetic  gases,  oxygen  is  paramagnetic  and  an  order  of  ten  to  a  thousand  times 
as  strong  as  the  diamagnetic  materials.  Oxygen  in  the  form  of  pressurized  air  will  be 
used  in  the  magnetic  fuel  cell  experiment  for  this  reason. 
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Figure  2  Gas  Pressure  versus  Susceptibility  for  Common  Gases  (Bismuth 

Included). 


(b)  Paramagnetism 


Air  is  the  predominant  paramagnetic  gas.  Comprised  of  19  percent 
oxygen  and  80  percent  nitrogen,  the  highly  paramagnetic  effect  j  =  1.87  x  10  ®)of 

oxygen  overcomes  the  small  diamagnetic  effect  ( j  =  -5.4  x  10  )  of  nitrogen. 
Paramagnetic  gases  are  temperature  dependant.  Hence,  heating  air  to  approximately 
300  °F  will  cause  oxygen  to  become  diamagnetic.  This  change  in  susceptibility  due  to 
temperature  or  pressure  will  be  key  in  achieving  the  optimal  fuel  cell  design.  The 
quantum  mechanical  equations  for  paramagnetism  are  given  in  Equation  (1.8)  [7,  8],  for 
two  situations.  The  first  situation  calculates  susceptibility  based  upon  changing  pressure 
with  a  fixed  temperature.  The  second  situation  calculates  susceptibility  based  upon 
increasing  temperature  [8].  Note  that  the  second  situation  shows  how  a  paramagnetic 
material  becomes  diamagnetic  upon  increased  temperature. 


7.1x10^'"  A 


K,T 


yPmH  j 


-csch 


V^bTj 


T  =  293  K 
T  >  293K 


(1.8) 


where; 


•  jU^  is  the  magnetic  moment  of  the  atom. 
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•  Kg  is  Boltzmann  constant  and 

•  r  is  temperature  (Kelvin). 


3,  Polarity 

Magnetie  fields  exhibit  a  permanent  dipole  state.  There  is  always  a  north  and 
south  pole.  Convention  states  that  the  north  pole  is  the  souree  and  the  south  pole  is  the 
sink.  Homogenous  magnetie  fields  are  a  key  ingredient  in  the  magnetie  fuel  eell.  Sinee 
diamagnetie  gases  will  eolleet  in  regions  of  lowest  magnetie  field  strength,  a  homogenous 
field  must  be  present  to  ensure  the  gas  stays  within  the  maximum  field  region. 


4,  Hysteresis 

For  anyone  wanting  to  work  with  magnetism,  understanding  magnetie  hysteresis 
is  essential. 


When  a  ferromagnetie  material  is  magnetized  in  one  direetion, 
(Figure  3)  it  will  not  relax  baek  to  zero  magnetization  when  the  imposed 
magnetizing  field  is  removed.  It  must  be  driven  baek  to  zero  by  a  field  in 
the  opposite  direetion.  If  an  alternating  magnetie  field  is  applied  to  the 
material,  its  magnetization  will  traee  out  a  loop  ealled  a  hysteresis  loop. 
The  laek  of  retraeeability  of  the  magnetization  eurve  is  the  property  ealled 
hysteresis  and  it  is  related  to  the  existenee  of  magnetie  domains  in  the 
material.  Onee  the  magnetie  domains  are  re-oriented,  it  takes  energy  to 
turn  them  baek  again.  This  property  of  ferromagnetie  materials  is  useful  as 
a  magnetie  "memory."  Some  eompositions  of  ferromagnetie  materials  will 
retain  an  imposed  magnetization  indefinitely  and  are  useful  as  "permanent 
magnets."  The  magnetie  memory  aspeets  of  iron  and  ehromium  oxides 
make  them  useful  in  audio  tape  reeording  and  for  the  magnetie  storage  of 
data  on  eomputer  disks.  [9] 
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Magnetization 


Material  magnetized 
to  saturation  by 
alignment  of  domains. 


Toward  saturation  in  dropped  to  zero  and  the  material  will 

the  opposite  direction  magnetization  (H 


remembers  its  history). 


Figure  3  Hysteresis  Curve  (from  [9]). 


5,  Magnetic  Circuits 

Magnetic  circuit  configurations  are  used  in  transformers  and  electrical  machines. 
Magnetic  circuit  analysis  techniques  are  analogous  to  those  of  electrical  circuits.  The 
magneto-motive  force  (mmf)  of  an  n-turn  current-carrying  coil  is  given  by  [5]: 


(1.9) 


where; 


3  is  the  magneto-motive  force, 


dL  is  change  in  unit  length, 


•  n  is  number  of  windings, 

•  /is  the  induced  current, 

•  is  core  magnetic  field  strength, 

•  is  core  length, 

•  is  gap  magnetic  field  strength  and 

•  is  gap  length. 

A  current-carrying  coil  is  the  magnetic  circuit  equivalent  to  a  voltage  source  in  an 


electrical  circuit.  Therefore,  magneto-motive  force  can  be  considered  to  be  analogous  to  a 
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source  voltage.  The  reluctance  (M)  of  a  path  for  magnetie  flux,  sueh  as  the  bar  or  iron  is 
given  by, 

91  =  ^,  (1.10) 

juA 

where  A  is  the  eross-sectional  area  of  the  magnetic  material.  Reluetance  is  analogous  to 
resistanee  in  an  eleetrical  eireuit.  When  the  bar  is  not  straight,  the  length  of  the  path  is 
taken  to  be  the  length  of  the  eenterline.  Therefore,  is  eonsidered  as  the  mean  length  of 
the  core  path.  Magnetic  flux  ( )  in  a  magnetic  circuit  is  equivalent  of  current  in  an 
electrieal  eireuit.  Magnetic  flux,  reluctance  ( 91 ),  and  magneto-motive  force  are  related 
by 

3  =  91^zi.  (1.11) 

This  is  analogous  to  Ohm’s  Law, 

V  =  IR,  (1.12) 

where; 

•  V  is  voltage  and 

•  i?  is  resistanee. 

The  advantage  of  the  magnetic  circuit  approach  is  that  it  can  be  applied  to 
unsymmetrieal  magnetic  cores  with  multiple  eoils.  Coils  are  sources  of  magnetomotive 
forces  that  can  be  manipulated,  in  the  same  way  voltage  sourees  are  in  electrieal  circuits. 
Reluctances  add  in  series  or  in  parallel,  just  the  same  as  resistance.  Fluxes  are  analogous 
to  eurrents.  The  magnetic  circuit  approach  is  not  an  exact  method  for  determining 
magnetie  fields.  However,  this  approach  is  sufficiently  precise  for  many  engineering 
applieations  [5]. 


B,  GAS  DYNAMICS 

Gases  at  high-pressure  travel  to  areas  of  low  pressure,  this  is  aeeording  to 
Bernoulli’s  prineiple  [1],  Equation  (1.13).  These  equations  ean  be  further  broken  down 
as  shown  in  Equations  (1.14)-(1.16), 
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where: 


^KE  +  ^PE  +  ^gas  ~  0  ’ 

(1.13) 

(1.14) 

AEpE  =  mg{h^  -  A.)  and 

(1.15) 

AE  =V(P,-P). 

gas  gV  /  • 

(1.16) 

•  ^ke  is  change  in  kinetic  energy  of  the  gas, 

•  ^PE  is  the  change  in  gravitational  potential  energy  of  the  gas, 

•  ^gas  is  the  change  in  pressure  energy  of  the  gas, 

•  is  the  final  velocity  of  the  gas, 

•  V.  is  the  initial  velocity  of  the  gas, 

•  g  is  gravity  (9.8  m/s^), 

•  hfis  final  height, 

•  m  is  the  mass  of  the  gas, 

•  /?.  is  initial  height, 

•  is  volume  of  the  gas, 

•  final  steady  state  pressure  of  the  gas,  and 

•  /’is  the  initial  pressure  of  the  gas. 

Re-arranging  Bernoulli’s  principle,  the  final  velocity  of  the  gas  can  be  solved  for 
as  seen  in  Equation  (1.17).  The  velocity  of  the  gas  is  proportional  to  the  square  root  of 
the  change  in  pressure  divided  by  the  density  change, 

llAP 

^1—’ 

where  v  is  gas  velocity  and  AP  is  change  in  pressure  of  the  gas. 
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(1.17) 


For  an  ideal  gas,  the  equation  of  state  [1]  is  Equation  (1.18).  Pressure  can  be 
determined  using  Equation  (1.19).  The  mass  of  the  gas  can  be  solved  using  Equation 
(1.20).  Einally,  using  Equations  (1.21)  and  (1.22),  the  gas  density  is  solved  for  in 
Equation  (1.22), 


PV,=n^RT, 

(1.18) 

n  RT 

(1.19) 

p  —  m 

^  Tr  ’ 

V 

m  =  n^M  , 

(1.20) 

II 

(1.21) 

PM 

(1.22) 

RT  ’ 

where:  is  the  number  of  moles  and  R  is  the  universal  gas  constant 

=  8.3 14  J/ (mol-K)) .  Since  the  gas  is  pulsated  into  a  chamber  at  close  to  one- 

millisecond  time  intervals  there  is  not  enough  time  for  heat  transfer  to  occur.  Therefore, 
this  process  can  be  analyzed  adiabatically  by  the  following  Equation  [1], 

EEJ  =  constant,  (1-23) 

where  y  is  the  ratio  of  molar  specific  heats  (y  =  1.4  for  diatomic  molecules). 


C.  FARADAY’S  LAW 

Earaday’s  law  states  that  the  line  integral  of  the  electric  field  around  a  stationary 
loop  equals  the  negative  surface  integral  of  the  time  rate  of  change  of  the  magnetic  flux 
density  B  integrated  over  the  loop  area  [10].  Eor  a  fuel  cell,  the  voltage  induced  reduces 
to 


^(|)  ^BA 

V  =  -n  =  -n - , 

At  At 


(1.24) 


where: 


A^zi  is  the  change  in  magnetic  flux, 
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•  is  the  change  in  time  and 

•  AS  is  the  change  in  magnetic  flux  density. 


D,  BASIC  CONCEPT 

From  the  previous  theory,  the  basic  equation  development  of  the  fuel  cell.  Figure 
4  shows  the  conceptual  operation  of  a  fuel  cell  is  presented.  (Note,  for  explanatory 
purposes  flux  lines  are  shown  doubling  in  the  soft  magnetic  core  vice;  in  reality,  the 
doubling  occurs  in  the  magnetic  fringe  around  the  magnets.)  At  time  zero,  the 
diamagnetic  gas  volume  is  outside  the  gap.  At  time  one,  the  diamagnetic  gas  enters  the 
gap  and  suppresses  the  local  magnetic  field.  Since  magnetic  flux  is  conserved,  the 
magnetic  flux  in  the  soft  magnetic  core  doubles.  This  change  in  flux  induces  a  voltage  in 
the  coil.  At  time  two,  the  diamagnetic  gas  volume  has  reduced  in  the  gap.  The  flux  is  no 
longer  sup-pressed,  thus  causing  the  flux  in  the  core  to  decrease  by  one  half  This  change 
in  flux  induces  a  negative  voltage  in  the  coil.  The  process  then  repeats  at  time  one. 


Soft  Magnetic 
Material  Core 

Magnetic 
flux 
in  gap 

N^S 

Magnet 


Voltage  Output 
Vout  =  0 


Diamagnetic 
Gas  Volume 


Voltage  Output 
Vout  >  0 


t=0 


t=l 


Voltage  Output 
Vout  <  0 


t=2 


•Cycle  starts 
•Magnetic  Flux 
doubles  in  core. 

•Change  in  flux  with  a 
change  in  time  produces 
a  voltage. 


•Magnetic  Flux 
halves  in  core. 

•Change  in  flux  with  a 
change  in  time  produces 
an  opposite  voltage. 
•Cycle  repeats  to  t=l. 


Figure  4  Conceptual  Operation  of  the  Magnetic  Fuel  Cell. 


Voltage  is  induced  by  a  change  in  flux  inside  a  coil  of  n  turns.  Equation  (1.25). 
A  change  in  flux  is  proportional  to  a  change  in  magnetic  flux  density  times  the  area. 


16 


(1.25) 


M  ^BA 

V  =  -n  =  -n - 

M  At 

The  magnetic  flux  density  is  given  by 

B  =  ^S\  +  Z)H.  (1.26) 

Therefore,  a  change  in  magnetic  flux  density  will  be  the  difference  between  the 
initial  susceptibilities  and  magnetic  field  strengths  and  final  susceptibilities  and  magnetic 
field  strengths, 

A5  =  //„(1  +  Zf)H,  -//„(!  +  z^)H,  «  B^HAz  ,  (1 .27) 

where: 

•  Zf  is  the  final  susceptibility, 

•  Zi  is  the  initial  susceptibility, 

•  Az  is  the  change  in  susceptibility, 

•  H is  the  final  magnetic  field  strength,  and 

•  H^  is  the  initial  magnetic  field  strength. 


For  this  basic  derivation  the  final  and  initial  magnetic  field  strengths  are  assumed  to  be 
equal.  This  will  not  be  assumed  for  the  advanced  derivation  shown  in  Chapter  111. 

Substituting  Equation  (1.27)  into  Equation  (1.25),  the  voltage  is  now  shown  to  be 
proportional  to  the  change  in  susceptibility, 

(1.: 

At 


The  unit-less  susceptibility,  z ,  is  shown  to  be  proportional  to  the  experimental 
susceptibility,  k  ,  by 


M 


Anp 


AnpK 

M 


(1.29) 


The  difference  in  susceptibility  using  Equation  (1.29)  becomes  Equation  (1.30)  and  is 
now  proportional  to  the  change  in  density  of  the  gas. 


17 


(1.30) 


Anp^K  Anp^K 
M  M 


Substituting  Equation  (1.30)  into  (1.28)  leads  to 


V  =  -n- 


At 


(1.31) 


where  voltage  is  proportional  to  the  ehange  in  density  of  the  gas.  The  density  of  the  gas 
ean  be  related  to  pressure 


PM 
RT  ■ 


(1.32) 


(The  assumption  that  the  gas  is  ideal  is  for  the  simplified  ease  and  is  treated  more 
rigorously  in  Chapter  111.)  Substituting  (1.32)  into  (1.31)  shows  voltage  is  proportional  to 
ehange  in  pressure  of  the  gas, 


V  =  -Ann 


^A  — 
RT  At  ■ 


(1.33) 


The  time  it  takes  the  gas  to  ehange  pressure  is  equal  to  the  height  of  the  fuel  eell  divided 
by  the  veloeity  of  the  gas  ( v), 

At  =  -,  (1.34) 

V 


where/?  is  the  height  of  the  gap. 

From  Bernoulli’s  prineiple  the  veloeity  of  gas  is  proportional  to  the  square  root  of 
two  times  the  ehange  in  pressure  of  the  gas  divided  by  the  density  of  the  gas.  Using 
Equation  (1.22),  density  ean  be  ehanged  into  equivalent  pressure,/’  to  give 


V  = 


2RT  AP 
M  Z’ 


(1.35) 


(Note:  Equation  (1.35)  has  been  simplified  for  the  basie  derivation;  a  more  eomplete 
analysis  is  in  Chapter  111.)  Substituting  Equations  (1.34)  and  (1.35)  into  (1.33)  produees 


V  =  -Ann 


(1.36) 
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Assuming  the  initial  pressure  is  mueh  greater  than  final  pressure  and  substituting  in  for 
area  ( ^4 )  of  the  fuel  eell  we  find 


V  =  -47rnuHwApJ^—  .  (1.37) 

V  MET 

As  seen  from  Equation  (1.37),  voltage  is  produeed  from  a  ehange  in  pressure  of  a 
gas.  Voltage  is  primarily  optimized  by:  1)  maximizing  the  pressure  ehange,  AP ;  2) 
maximizing  the  number  of  turns  of  the  eoil,  n  ;  3)  maximizing  width  of  the  fuel  eell,  w ; 
and  4)  maximizing  the  magnetie  field  strength  in  the  gap,  H . 

This  Chapter  has  introduced  the  basic  concepts  required  to  understand  how  a 
magnetic  fuel  cell  functions.  Extensive  discussion  of  magnetic  properties  was  reviewed, 
along  with  gas  dynamics.  Using  these  concepts,  the  basic  equation  for  the  fuel  cell  was 
derived  showing  its  dependency  on  gas  pressure,  magnetic  field  strength,  and  fuel  cell 
geometry.  Chapter  III  will  present,  a  more  thorough  analysis  of  the  magnetic  fuel  cell,  a 
discussion  and  derivation  of  coil  optimization,  and  what  the  implementation  issues  were 
for  the  experiment. 
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III.  EXPERIMENTAL  SETUP 


A,  THEORETICAL  DEVELOPMENT 

This  Chapter  will  discuss  the  detailed  derivation  of  the  magnetic  fuel  cell.  Two 
key  parameters  that  will  be  solved  are  the  static  magnetic  field  inside  the  pressure 
chamber  and  the  change  in  magnetic  field  due  to  the  pressure  change.  These  parameters 
affect  the  magnitude  of  the  voltage  in  the  fuel  cell  equation.  The  next  Section  discusses 
the  optimization  of  power.  The  final  Section  of  Chapter  III  will  apply  the  detailed 
derivation  to  the  experiment  model. 

The  geometric  configuration  of  the  magnetic  fuel  cell  is  shown  in  Figure  5.  A 
stable  DC  magnetic  field  is  generated  using  Nd-Fe-B  permanent  magnets.  These 
magnets  are  placed  on  the  ends  of  highly  permeable,  high  frequency,  soft  magnetic 
material  (called  a  c-core).  A  pressure  chamber  fits  inside  the  gap  between  the  magnets. 
High-pressure  gas  flows  into  and  out  of  the  cylinder  as  controlled  by  valves  two  and 
three.  The  change  in  magnetization  of  the  gas  in  the  gap  is  caused  by  a  change  in 
pressure  of  the  gas.  The  change  in  magnetization  causes  a  change  in  magnetic  flux.  This 
change  in  magnetic  flux  induces  a  voltage  in  a  coil  wound  around  the  c-core  soft 
magnetic  material.  The  induced  voltage  is  applied  to  the  load  resistance. 


Magnetic  Fuel  Cell:  Gas  Flow  Path  Magnetic  Fuel  Cell:  Magnetic  Circuit 

Pressure 


Ps 

Figure  5  Magnetic  Fuel  Cell  Gas  and  Magnetic  Conceptual  Flow  Paths. 
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The  first  part  of  the  derivation  is  a  three-step  development  to  determine  the  time 
the  gas  takes  to  enter  the  ehamber  and  exit  the  ehamber.  (Note:  that  the  viseosity  effeets 
and  heating  effeets  of  the  gas  are  ignored  in  this  derivation  due  to  the  small  length  of 
piping  the  gas  travels  through  and  the  short  time  the  gas  takes  to  diseharge  to  the 
atmosphere.)  The  first  step  is  to  determine  the  time  for  the  gas  to  enter  the  gas  ehamber 
using  Bernoulli’s  prineiple, 

P  ^  p  ^ 

_ ^  ^ valve!  j 

Ps  2  Pgap'i  ^ 

where: 

•  ^vaivei  is  veloeity  of  the  gas  entering  valve  two, 

•  is  the  veloeity  of  the  gas  in  the  supply  tank, 

•  is  the  pressure  of  the  gas  in  the  supply  tank, 

•  is  the  density  of  the  gas  in  the  supply  line, 

•  is  the  pressure  of  the  gas  in  the  gap,  and 

•  Pgap  is  tiiG  density  of  the  gas  in  the  gap. 


Using  Bernoulli’s  prineiple  [1],  the  veloeity  of  the  gas  traveling  through  valve  two  is 
determined  by. 


V 


valve! 


(2.2) 


Next,  the  final  veloeity  of  the  gas  entering  the  pressure  ehamber  is  determined  using  the 
eontinuity  equation  [1], 

P s'^valvel^ valve!  P gap'^gap^ gap-in  ’ 


where: 


•  Aaive!  ’  eross-sectional  area  of  valve  two, 

•  ,  the  eross-seetional  area  of  the  gas  ehamber  and 
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•  ^gap-in  ’  the  velocity  of  the  gas  entering  the  gap. 


The  supply  volume,  ,  is  much  larger  than  the  gap  volume,  Vgap.  The  density  of  the  line 
equals  the  density  of  the  gap.  Rearranging  Equation  (2.4)  produces 


The  time  for  the  gas  to  enter  the  gas  chamber  is. 


= 


where  is  the  height  of  the  gas  chamber.  Substituting,  Equation  (2.5)  and  (2.2)  into 
(2.6)  produces 


At,  „  = 


A  P  P 

^valvel  2  ^ _ 

^gap  \  K  Ps  P gap 


The  second  step  is  to  determine  the  time  for  the  gas  to  exit  the  gas  chamber. 
Applying  Bernoulli’s  principle  and  solving  for  the  velocity  of  the  gas  entering  valve 
three. 


P  p 

gap  '^gap-out  _  a  gm  '^valvei 


P  P 

...  _  S^P  _  ^  atm 

^valve?,  ~  A  ^ 

V  V  ^gaP  P  atm 


where: 


•  ^vaiveT,  ^hc  vclocity  of  the  gas  entering  valve  three. 


•  ^gap-out  is  velocity  of  the  gas  leaving  the  gap. 


Palm  is  atmospheric  pressure,  and 


Pam  is  atmospheric  density. 
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Next,  the  final  velocity  of  the  gas  exiting  the  pressure  chamber  is  determined  using  the 
continuity  equation, 


P atm'^valve3^atm 


P 2ap'^^an^valve3  ’ 


(2.10) 


where  is  the  cross  sectional  area  of  valve  two  and  is  the  velocity  of  the  gas 
exiting  valve  two.  Rearranging  terms  and  solving  for  atmospheric  velocity. 


V  = 

atm 


PeaD^miP 


gap  gap*^  valve?) 


PatA 


atm'^  \alve? 


The  time  for  the  gas  to  exit  the  gas  chamber,  ,  is. 


Ku,=- 


Substituting,  Equation  (2.9)  and  (2.1 1)  into  (2.12)  gives 


Kur  = 


0  A 

h'  atm^valve? 


h 


gap 


P 


A  P  ^ 

gap  _  atm 
P gap  P atm  J 


(2.11) 


(2.12) 


(2.13) 


The  third  step  is  to  determine  the  total  time  of  operation  of  the  magnetic  fuel  cell. 
This  process  will  include  the  time  taken  to  cycle  open  and  shut  valves  two  and  three.  The 
cycle  time  is  the  sum  of  the  valve  operation  times  and  charge  and  discharge  times, 

^ cycle  ^ valve!  +  Kut  +  ^ valve?  ’  (2.14) 

where; 

•  ^ cycle  is  cycle  time  of  the  magnetic  fuel  cell, 

•  Kaive  is  tii^^  io  open  and  shut  valve  two,  and 

•  C/vesis  the  time  to  open  and  shut  valve  three. 


The  cycle  frequency,  ,  is 


f cycle 


^  cycle 


(2.15) 
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The  volume  flow  rate  of  the  eyele,  ,  is 


4w,=^.  (2.16) 

^  cycle 

The  duration  or  total  time  it  takes  for  the  supplying  high-pressure  gas  tank  to  beeome 
exhausted,  ,  is 


V, 


r 


cycle 


(2.17) 


The  timing  of  the  fuel  eell  has  been  determined.  The  next  series  of  steps  analyzes  the 
magnetie  eireuit  (Figure  5)  of  the  magnetie  fuel  eell.  This  analysis  begins  with  the 
determination  of  the  magnetie  load  line.  From  the  statie  magnetie  field  equations,  the 
eireulation  of  magnetie  field  strength  per  unit  length  is  equal  to  the  integral  of  the  eurrent 
density  per  unit  area, 

^HdL  =  ^^Jds  =  nI ,  (2.18) 


where  J  is  the  indueed  eurrent  density  in  the  set  of  windings,  andds  is  the  unit  area 
perpendieular  to  the  magnetie  flux  lines.  The  eireulation  of  the  magnetie  flux  density 
around  a  elosed  surfaee  is  zero, 

^B-ds=d.  (2.19) 

5 


(Assumption;  the  flux  generated  by  the  field  souree  is  eonfined  to  the  magnetie  eireuit 
with  minimal  leakage.)  Equation  (2.20)  beeomes 

+  ^gap^gap  +  ^ coreLcore  =  ’  (2.20) 

where: 


•  is  the  magnet  magnetie  field  strength. 


•  is  the  length  of  the  magnets. 


•  ^gap  is  the  gap  magnetie  field  strength. 
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•  is  the  length  of  the  gas  chamber  or  gap, 

•  ^core  is  the  core  magnetic  field  strength,  and 

•  is  the  length  of  the  core. 

The  magnetic  flux  through  a  surface  is  zero  as  shown  by  Equation  (2.19).  The  magnetic 
flux  densities  for  each  component  are  therefore  equal  (assuming  fringing  flux  is  small), 

B  =  B  =B  =B  ,  (2.21) 

gap  core  m  ^  V  / 

where; 

•  B^^^  is  the  gap  magnetic  flux  density, 

•  is  the  core  magnetic  flux  density  and 

•  B^  is  the  magnet  magnetic  flux  density. 

The  basic  equation  for  magnetic  field  strength  becomes  (Note:  the  derivation  of  Equation 
(2.25)  is  covered  in  Appendix  B.), 


H  = 


B 


B 


TT  _ 

core 


B 


B core  Bo 


Bo 


k  = 


nwh„ 


2fi  Jh^  + 

t^core  \  c  c 


where; 


•  Bgap  is  the  gas  permeability  within  the  gas  chamber. 


(2.22) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 
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•  l^core  is  permeability, 

•  A:  is  a  simplification  constant, 

•  is  the  core  height,  and 

•  is  the  core  width. 

Assuming  the  magnet  has  a  linear  second  quadrant  demagnetization  curve,  then  the 
equation  is  of  the  form  [5], 

{221) 

and  the  magnetic  permeability  is  defined  by, 

(2.28) 

where; 

•  is  the  residual  magnetization, 

•  /u^  is  the  permeability  of  magnet  and 


•  is  the  magnet  coercivity. 


Setting  magnetic  flux  densities  equal  to  each  other.  Equations  (2.27)  and  (2.28)  allow  the 
calculation  of  the  magnetic  load  line.  Substituting  for  magnetic  field  strength  from 
Equations,  (2.22)-(2.25),  into  Equation  (2.20)  produces  Equation  (2.29), 


-HL 


2B  ^  B  ^  Bk 

- L  H - L  -n  — 

P 

MgapMc  PcorePo  Po 


(2.29) 


Rearranging  Equation  (2.29)  to  solve  for  magnetic  flux  density  produces. 


B  =  - 


L 


gap 


L 


(2.30) 


-nk 


P 


gap 


Pc, 
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Setting  Equation  (2.27)  equal  to  Equation  (2.30)  finds  the  magnetie  field  strength  of  the 
magnets, 


D  I 


= 


f^gap 

-B.. 


+  ^^-nk 

Bcore 


BoK 


+ 


L 


+ 


-nk 


A 


Bg 


Me. 


The  flux  through  a  surface,  (j) ,  is  given  by 


(2.31) 


(/)  =  ■  ds  . 

s 


(2.32) 


The  area  of  the  circuit  is  fixed.  Equation  (2.32)  will  reduce  to 

(l)  =  BA,  (2.33) 


where  A  is  the  cross-sectional  area  of  magnets. 

E araday’s  Induction  equation  describes  the  voltage  induced  into  the  windings  by 
the  change  in  flux  in  the  gas  chamber  (repeated  here  for  convenience). 


V  =  -n 


At 


(2.34) 


Substituting  magnetic  flux  density  from  the  load  line  equation  into  Earaday’s  Equation 
produces  Equation  (2.35)  for  the  gas  entering  the  chamber  and  Equation  (2.36)  for  the  gas 
discharging  from  the  chamber. 


B,=- 


MoH„,K 


L 


-  +  - 


L 


gap  aore__^J^ 


Mgapl  Me, 


D  _ _ 

out  -  I 


MoH^K 


L 


gap 

Mgapl  M, 


+  ^^-nk 


(2.35) 


(2.36) 
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where  B.^  is  the  magnetic  flux  density  at  high  pressure  and  is  the  magnetic  flux 
density  at  atmospheric  pressure. 


The  change  in  magnetic  flux  density  becomes, 


AB  =  - 


L 


-  +  - 


L 


■  +  ■ 


uH  L 

r^o  m  m 


-nk 


L 


gap 


+ 


L 


(2.37) 


-nk 


Mgap2  Bcore  Bgapl  Bcore 

Substituting  Equation  (2.37)  into  (2.34)  and  using  (2.33),  the  voltage  can  be  determined. 


V  = 


1 

At, 


ly  Bgapl  Bcore  Bgapl  Bcore 


(2.38) 


Using  Equations  (2.7)  and  (2.13)  to  substitute  for  time  produces 


V  =  PjrH,^A 


\alve2 


gap 


'  gap 


Ps  P. 


PgapA 


gap  J 


gap  gap 


P alm^valve2  ^ 


gap 

\Pgap 


P, 


atm  J 


L 


gap 


+ 


L 


-nk 


gap 


+  - 


L 


-nk 


Pgap\  Pa 


Pgapl  Pc. 


.  (2.39) 


The  following  simplification  will  help  to  reduce  the  complexity  of  this  equation 
producing  Equation  (2.40): 


P  p 

•  ^^  =  ^^forA4, 

P gap  P atm 

•  ^  =  —  for  ^Kut  and  V  «  , 

Pgap  Ps 

•  Patm  *  1 ,  and 


•  ^  =  K^g- 
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To  allow  a  full  theoretical  analysis  the  following  substitution,  =  J2 
made. 


s  atm 


V  Ps  Patm  ) 


IS 


( 


V  =  unHLwy 

r^o  mm2. 


g  gas 


\alve2 


A 


gap 


I 

'valve'i 


s  gap 

4. 


^Qnn  L. 


gap  I  core 


-k 


L. 


gap  _j_  core  _ 


^gap\  t^core  ^  gapl  f^core  ) 

The  electrical  energy  output  of  the  magnetic  fuel  cell  is 

_ 

cycle  ~  ^ cycle  ’ 


(2.40) 


(2.41) 


where  ^^yde  electrical  energy  produced  by  the  magnetic  fuel  cell  each  cycle,  and  R  is 
load  resistance.  The  power  of  the  magnetic  fuel  cell  is 

Pceii=—,  (2.42) 

where  is  the  power  produced  by  the  magnetic  fuel  cell  each  cycle. 

To  determine  efficiency,  first  an  energy  balance  will  be  used  to  balance  the 
energy  removed  from  the  gas  to  the  energy  produced  by  the  fuel  cell  (assuming  no 
losses), 

K  AP  =AE  ,  .  (2.43) 

gap  mag  cycle  V  ' 

This  becomes 


V^t 

AD  _  ayala 
mag  n  y  . 

P^gap 


(2.44) 


where  is  the  magnetic  pressure  drop  caused  by  the  magnetic  fuel  cell,  and  is 

the  volume  of  the  gap.  The  efficiency  is  derived  by  using  the  energy  of  the  fuel  cell 
divided  by  the  total  energy  of  the  gas. 


s  = 


AP 


cell 


AP, 


total 


(2.45) 
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(2.46) 


_ r!i _ 

PV -P  (V +V  )’ 

where; 

•  s ,  efficieney  of  the  magnetie  fuel  eell, 

•  total  ’  total  energy  of  the  gas,  and 

•  l^tteii  ’  total  energy  of  the  magnetic  fuel  cell. 

The  efficiency  equation  can  be  simplffied  by  assuming, 

PV  »  P  (V  +V  ) 

s  atm\^  s  ^  ^  gap)  ’ 


therefore  , 


s  = 


V% 
RPsK 


Several  substitutions  can  be  made  shown  below: 
•  substitute  in  duration  time. 


RPsKyole  ’ 


(2.47) 


(2.48) 


(2.49) 


•  substitute  in  volume  flow  rate. 


V\ 


s  - 


cycle 


RPsV,a, 


•  substituting  into  Equation  (2.50),  the  effect  of  magnetic  pressure, 

1 


V^t 


E  = 


cycle 


V^t 


cycle 


R{P-^P  \V 

\  s  ^magj'g 


f 


gap 


R 


P  - 


V^t. 


cycle 


RV 


gap  J 


(  RP.Kao 

F 

s  gap 

1 

gap 

y 

^  ^  cycle 

y 

(2.50) 


(2.51) 


•  substitute  in  voltage. 
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s  = 


^  cycle 


RPV 

s'  gap 


LI  nH  L  w  V 

f^o  m  m  g  gas 


\alve2 


Ps\ 

A 


'valve2 


-k 


^(Ton 


-k 


yPgapl  Pc, 


Pgapl  Pc, 


,(2.52) 


JJ 


substitute  in  volume  of  gap, 


^Sap  ^gap^gap^gap  ’ 


(2.53) 


£  = 


W  t  , 

gap  cycle 


RPh  L 

s  gap  gap 


Po^H^L^v 


m  m  gas 


A. 


Kalvel 


A 


gap 


^gap  L, 


pA 

A 


’give  3 


'core 


L 


T. 


W 


.(2.54) 


JJ 


From  Equation  (2.54),  several  observations  ean  be  made.  The  magnetic 
permeability  of  free  space,  /u^ ,  times  the  magnetic  field  strength,  of  the  magnet  is 
approximately  one.  The  number  of  turns  variable,  n  ,  is  squared  but  there  is  an  implicit 
relationship  between  n  and  load  resistance,  R  .  The  impact  of  this  is  that  maximizing 
n  will  also  raise  resistance  R  when  trying  to  maximize  power  transfer.  This  will  be 
discussed  in  detail  in  Section  B.  The  magnetic  terms  under  the  denominator  in  Equation 
(2.54),  -  A: ,  are  the  same  for  both  denominators  with  the  exception  of  the 

permeability  of  the  gap  terms,  and  p^^^^  ■  Therefore,  to  maximize  voltage  and 
efficiency  both  the  core  length  divided  by  the  core  permeability,  /  p^^A  ~  ^ 

k  term  should  be  minimized  using  a  large  core  permeability  value,  p^^^^ .  The  velocity  of 
the  gas,  also  has  a  dependent  relationship  between  itself  and  supply  tank  pressure, 

;  the  consequence  being,  minimizing  supply  tank  pressure  also  minimizes  gas  velocity. 

The  rest  of  the  terms  reveal  order  of  magnitude  effects  on  efficiency  based  upon  their 
relative  magnitude  and  are  listed  in  Equations,(2.55)-(2.59),  below,  in  order  of  decreasing 
importance: 

^oc^,  (2.55) 

gap 

£°cLl,  (2.56) 
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s  ^w, 


(2.57) 


s 


(2.58) 

(2.59) 


B,  RESISTANCE  ESTIMATION  AND  COIL  OPTIMIZATION 

The  number  of  windings  term,  n  ,  in  Faraday’s  Induetion  Equation,  (2.34),  is  a 
key  parameter  effecting  coil  optimization.  On  the  surface,  a  larger  number  of  windings 
should  generate  a  large  voltage  output  (Equation  (2.34)).  But,  looking  at  power  produced 
by  the  fuel  cell  shows  that  internal  resistance  and  therefore  load  resistance  also  increases 
as  number  of  windings  increases  (Equation(2.60)).  Increasing  the  number  of  windings 
causes  the  internal  current  to  drop  and  therefore  output  power  to  drop.  This  Section 
optimizes  the  power  factor  term,  /7? ,  in  the  power  equation  with  respect  to  winding 
resistance,  since 


V  =  -n - 

At 


R 


P  cc 


(2.60) 


In  order  to  establish  equivalency  between  two  different  gauge  size  wires,  a  fixed 
winding  cross-section  is  used  (refer  to  Eigure  6)  with  the  horizontal  distance  represented 
as  X  and  the  vertical  distance  represented  as  y  .  This  analysis  is  based  upon  two 
simplifying  assumptions,  one  being  the  wire  is  wound  around  a  circular  magnetic  core 
with  radius,  a  ,  and  the  other  being  the  wire  itself  is  circular  with  radius,  b  . 


Figure  6  Two  Examples  of  Wire  Winding  Cross-Sections. 


The  number  of  windings  fitting  onto  each  layer,  ,  of  the  coil  can  be  determined  using 
geometry. 
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^layer 


wire  diameter  2b 


(2.61) 


The  total  number  of  layers,  n  ,  is  then  determined  as 


"  2b 


(2.62) 


The  total  number  of  windings,  n  ,  is 


n  =  n,„n„  = 


xy 


layer  y  ^^2 


(2.63) 


The  length  of  the  first  layer,  ,  is  determined  using 

kayerx  =  (looP  circumference)  =  2nan,^^^^  (2.64) 

The  length  of  the  seeond  layer,  2  >  is  found  as 

Llayer  2  =  +  2b)  .  (2.65) 

The  total  length  of  the  coil,  ,  accounting  for  all  layers  is 

k=n^. 

kail  =  27:ni^yer  Z  (^  +  ^^(^  - 1))  .  (2.66) 

k=\ 

Substituting  in  Equations  (2.61)  and  (2.62)  in  Equation  (2.66)  produces 

k=^ 

X 

(2.67) 

^  k=\ 

The  equation  for  resistance  of  a  wire  [10]  is  given  by 

R  =  PjycoiL^  (2.68) 

'^wire 

where; 

•  is  the  resistivity  of  the  wire, 

•  L^^-i  is  the  total  length  of  the  coil,  and 
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•  is  the  cross  sectional  area  of  the  wire. 

Substituting  Equation  (2.67)  into  Equation  (2.68)  produces 

k=-^ 

R  =  ^J^{a  +  2b{k-\)).  (2.69) 

Tib 


Einally,  substituting,  Equations  (2.63)  and  (2.69)  back  into  the  power  factor  constant, 
Equation  (2.60),  produces  (normalizing  x,y  =  \) 


n 

R 


k= — 
2h 


16yC>b^  (a  +2b{k-V)) 


k=i 


(2.70) 


Equation  (2.70)  is  plotted  in  Eigure  7  with  a  varied  wire  cross-section,  b  ,  and  a  fixed 
loop  radius,  a  =  0. 1025  m .  This  plot  shows  that  changes  in  the  wire  cross-section  are 
almost  linearly  compensated  for  by  the  summation.  Therefore,  to  maximize  the  power 
factor  constant  in  Equation  (2.70),  minimize  the  loop  radius,  a  .  Equation  (2.70)  can  be 
simplified  to 


R  a 


(2.71) 
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Figure  7  Plot  of  Power  factor  Constant  for  Varying  Wire  Cross-Sections. 
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Using  actual  wire  gauge  sizes,  a  plot  of  power  versus  number  of  windings  is 
shown  in  Figure  8.  Faraday’s  Induction  law,  Equation  (2.72),  was  used  to  build  Figure  8. 
The  change  in  magnetic  flux  was  10^^  T  .  The  area  was  0.02  m^ .  The  change  in  time 
was  10^^  s  .  The  resulting  equation  was  then  found  to  be. 


J/2  _  1 

f  n^BA^ 

2  2 

n 

( 1  X  0.02  ^ 

400n" 

1 

K  At  j 

R 

1 

o 

R  ’ 

(2.72) 


where  7?  is  a  function  of  n.  Figure  8  shows  there  are  an  optimal  number  of  windings  to 
achieve  maximum  power.  This  plot  also  shows  that  lower  gauge  wires  produce  a  larger 
power  output  even  though  they  have  less  windings.  The  theory  for  the  magnetic  fuel  cell 
is  now  applied  to  the  experimental  setup. 


18  Gauge 
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0  200  400  600  800  1000 

Number  of  Windings 

Figure  8  Power  vs.  Number  of  Windings  for  Different  Gauge  Wires. 

C.  IMPLEMENTATION 

1,  Configuration  Decision 

The  equipment  determination  for  the  experiment  was  based  upon  low-cost 
commercial  components.  The  pressure  of  2000  psi  used  in  the  experiment  was  based 
upon  the  smallest  high-pressure  sample  chamber  offered  by  Sunnyvale  Valve  and  Fitting 
Company  [11].  The  smallest  chamber  was  also  chosen  to  minimize  the  length  of  the  gap. 
The  smallest  chamber  diameter  was  two  inches.  The  quarter-inch  copper  tubing  was 
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chosen  based  upon  being  already  available  on  site.  This  dictated  the  sample  chamber  and 
valve  inner  diameter  sizes  as  one-quarter  inch. 

The  permanent  magnets  were  obtained  from  Stanford  Magneties  Corporation 
[12].  Nd-B-Fe-40  magnets  having  a  very  high  residual  magnetic  field  and  high  eoereivity 
are  sold  at  a  low  cost.  The  number  and  size  of  the  magnets  were  chosen  based  upon  the 
height  and  diameter  of  the  sample  ehamber  (i.e.,  eight  inches  by  two  inches).  A  holding 
apparatus  for  the  magnets  was  designed  and  fabricated  by  George  Jaksha,  the  NFS 
physies  department  maehinist.  This  apparatus  held  the  permenant  magnets  against  the 
sample  chamber  and  maintained  alignment  with  respeet  to  eaeh  other.  An  air-gap 
magnetie  return  path  was  ehosen  over  a  soft  magnetie  material  core  path  due  to  excessive 
costs  of  the  soft  magnetie  materials. 

Due  to  the  high-speed  high-pressure  requirements  of  the  valves,  Swagelock  ball 
valves  were  chosen  to  provide  quick  operation  in  both  opening  and  closing.  Twenty-four 
gauge  wire  was  ehosen  for  the  magnetic  fuel  cell  coil  due  to  it  being  available  on  site  and 
in  large  quantity.  This  wire  allowed  a  large  number  of  windings  to  be  wrapped  around 
the  sample  ehamber  support  apparatus. 

2.  Experiment  Predictions 

With  the  material  geometries  established,  the  predicted  quantities  of  the 
experiment  ean  be  determined.  This  experiment,  being  the  first  of  its  kind,  foeused  on 
obtaining  an  output  voltage  from  the  fuel  cell.  Only  the  diseharge  portion  of  the  fuel  cell 
eycle  is  developed  sinee  this  theoretically  should  produce  the  largest  voltage.  The 
derivation  is  shown  in  Appendix  C.  The  input  values  and  constants  are  summarized  in 
Table  2  of  Appendix  D.  All  the  derived  values  are  shown  in  Table  3  of  Appendix  D.  For 
this  configuration,  the  predicted  voltage  to  be  sensed  was  -2.75  millivolts  to  occur  in 
roughly  3.66  milliseconds. 

This  chapter  performed  an  indepth  analysis  of  the  theory  of  the  magnetic  fuel  eell. 
A  brief  discussion  was  provided  on  optimizing  the  number  of  windings  of  the  magnetic 
fuel  cell  eoil.  Several  experimental  design  deeisions  were  diseussed  along  with  the 
predicted  voltage  response  of  the  experiment.  Having  determined  the  theoretical 
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characteristics  of  the  magnetic  fuel  cell,  Chapter  IV  discusses  the  experimental  setup  and 
implementation  of  the  experiment. 
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IV.  EXPERIMENT 


A,  OVERVIEW 

The  experiment  was  eonducted  in  the  Rail  Gun  Lab  in  Spanagel  Hall  at  the  Naval 
Postgraduate  Sehool  in  November  of  2002.  The  goal  of  the  experiment  was  to  verify 
theoretical  equations  and  predictions  developed  in  Chapter  III. 

The  magnetic  fuel  cell  is  shown  in  Figure  9.  The  picture  on  the  left  is  without  the 
coil  and  the  picture  on  the  right  is  with  the  coil.  In  the  picture  without  the  coil,  the  gray 
structure  wrapped  with  tygon  tubing  is  the  support  apparatus  mentioned  in  Chapter  III. 
Four  wood  inserts  were  provided  for  additional  support.  The  gas  tubing  with  connectors 
can  be  seen  on  either  end  of  the  pressure  chamber,  which  is  located  inside  the  support 
apparatus.  The  picture  on  the  right  of  Figure  9  is  the  fuel  cell  with  515  windings  of  25- 
gauge  wire  wrapped  around  it.  Duct  tape  was  used  to  ensure  the  windings  stayed  tightly 
wound  around  the  support  apparatus. 


Figure  9  Magnetic  Fuel  Cell:  (left)  without  coil,  (right)  with  coil. 


B,  SETUP 

A  diagram  of  the  experimental  setup  is  shown  in  Figure  10.  Nitrogen  and  air 
were  the  diamagnetic  and  paramagnetic  gases  used.  An  oscilloscope  was  used  to 
determine  output  voltage.  An  open  microphone  was  used  to  provide  an  audio  trigger 
caused  by  the  vibration  noise.  There  was  considerable  vibration  noise  developed  by  the 
magnetic  fuel  cell  due  to  the  rapid  expansion  of  2000-psi  gas  to  atmospheric  pressure. 
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Eight  two-inch  by  two-inch  by  one-quarter-inch  Nd-B-Fe-40  magnets  were  used 
in  the  experiment  to  magnetize  the  paramagnetie  gas.  The  magnets  eharacteristies  are 
listed  in  Table  1 .  An  explanation  of  eaeh  term  is  provided  below. 

•  The  residual  induetanee,  By,  is  the  value  of  induction  at  the  point  of  the 
hysteresis  loop,  at  whieh  the  hysteresis  loop  erosses  the  B  axis  at  zero 
magnetizing  foree.  The  By  represents  the  maximum  magnetie  flux  density 
output  of  the  magnet  without  an  external  magnetie  field. 

•  The  eoereivity  of  the  magnet,  He,  is  the  demagnetizing  foree,  measured  in 
Oersted,  neeessary  to  reduee  observed  induetion,  B  to  zero  after  the  magnet 
has  previously  been  brought  to  saturation. 

•  The  intrinsie  eoereivity,  Hd,  is  an  intrinsie  ability  of  a  magnet  to  resist 
demagnetization.  Its  value  is  measured  in  Oersted  and  eorresponds  to  zero 
intrinsie  induetion  in  the  magnet  after  saturation. 

•  The  (BH)max  term  whieh  is  a  representation  of  the  overall  energy  of  the 
magnet.  There  is  a  point  on  the  hysteresis  loop  at  whieh  the  produet  of 
magnetizing  foree  H  and  induetion  B  reaehes  a  maximum.  The  maximum 
value  is  ealled  the  Maximum  Energy  Product.  At  this  point,  the  volume  of 
magnet  material  required  to  projeet  a  given  energy  into  its  surrounding  is  a 
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minimum.  This  parameter  is  generally  used  to  describe  the  “strength"  of  the 
permanent  magnet  material. 

•  The  Curie  temperature  is  the  temperature  at  which  the  parallel  alignment  of 
elementary  magnetic  fields  completely  disappears  and  the  material  is  no 
longer  able  to  hold  magnetization. 

•  The  maximum  operating  temperature  is  the  temperature  above  which  the 
magnets  will  start  demagnetizing. 

Table  1.  Magnetic  Characteristics  ofNd-B-Fe-40  (after  [12]). _ 


Grade 

Br 

(KG) 

He 

(KOe) 

H,i 

(KOe) 

(BH)^ax 

(MGOe) 

Curie 

Temp. 

(“C) 

Max.  Op. 
Temp. 
(“C) 

N40 

12.5-12.8 

>11.6 

>12 

38-41 

310 

80 

Five  hundred  24-gauge  copper  windings  were  wound  around  the  pressure  chamber  to 
detect  the  magnetic  field  change  caused  by  the  pressurized  gas.  The  gas  cylinder 
containing  2000  psi  was  used  as  the  gas  source  for  the  experiment.  A  5600B  series  high- 
pressure  regulator  regulated  downstream  gas  pressure  to  1800  psi.  All  three  valves  used 
in  the  system  were  high-pressure  rated,  one-quarter-inch,  Swagelock  “40”  series  ball 
valves  [11]. 

Valve  one  acted  as  a  two-valve  isolation  between  the  high-pressure  tank  and 
downstream  piping.  Valve  two  controlled  gas  flow  into  the  cell  gap  and  isolated  the  gap 
from  upstream  pressure  when  valve  three  was  opened.  Valve  three  isolated  the  gap  from 
atmosphere  when  valve  two  was  open.  (Ideally,  an  automated  system  using  solenoid- 
operated  valves  would  charge  and  discharge  the  cell  gap  at  a  user-controlled  frequency.) 

The  pressurized  chamber  containing  a  two-inch  cell  gap  is  a  Swagelock  316L 
stainless  steel  sample  cylinder  with  an  eight-inch  length  and  a  one-quarter  inch  cell  wall 
rated  for  5000  psi  [11].  Both  valves  and  sample  chamber  were  chosen  to  be  stainless 
steel  316L  due  to  this  material’s  anti-ferromagnetic  properties  and  ability  to  handle  high- 
pressures.  The  piping  between  the  regulator  and  valve  one  is  one-quarter  inch  copper 
tubing.  The  piping  between  valves  one  and  two  was  7000  psi  rated  flexible  tubing. 

Valve  three  vented  the  high-pressure  gas  in  the  gap  to  atmosphere. 

The  eight  magnets  are  held  in  place  against  the  sample  chamber  by  a  fabricated 


PVC  cubic  container  with  rounded  wood  wedges  in  the  corners  (Figure  11).  Figure  1 1 
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shows  a  cutout  cross-section  of  the  magnetic  fuel  cell  with  dimensions.  A  Radio  Shack 
omni-microphone  was  used  to  trigger  the  scope  on  discharge  of  the  gas.  The  windings 
and  the  oscilloscope  were  connected  by  a  coaxial  cable. 


Wood 

Wedges 

0.24inch 


Figure  11  Cross  Section  of  Fuel  Cell. 

2,  Calibration 

A  Hall-effect  Gauss-meter  was  used  to  determine  the  magnetic  field  within  the 
de-pressurized  cylinder.  At  the  center-point  of  the  cylinder  with  the  probe  extending  in 
four  inches,  the  magnetic  field  was  measured  to  be  0.1  Tesla.  This  agreed  with  the 
predicted  value  of  0.1  Tesla  derived  in  Appendix  D.  To  verily  operation  of  the  outer  coil 
and  determine  the  coupling  parameter  of  the  circuit,  the  following  PSPICE  test  circuit 
was  simulated  (Figure  12).  A  seven-coil  wire  loop  was  made  from  24-gauge  insulated 
copper  wire  and  inserted  inside  the  pressure  chamber.  A  high  current  voltage  source  and 
high  power  IGBT  (Figure  13)  was  used  in  parallel  with  a  square  wave  generator  and  a 
voltage  source.  The  IGBT  was  triggered  off  a  square  wave  pulse  source  with  a  pulse 
width  of  30  milliseconds.  The  rise  time  on  the  pulse  was  nine  microseconds.  A  current 
pulse  of  two  amps  was  applied  to  the  loop.  The  output  of  both  the  pulse  source  and  fuel 
cell  windings  was  connected  to  channels  one  and  two  of  the  oscilloscope. 
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Figure  13  Voltage  Source  and  High  Power  IGBT. 


Both  the  test  signal  and  PSPICE  model  output  are  shown  in  Figure  14  and  Figure  15.  A 
significant  effort  was  made  to  get  the  test  signal  output  of  the  oscilloscope  and  the 
PSPICE  model  output  to  match.  The  test  coil  and  fuel  cell  coil  were  measured  for 
inductance.  The  test  signal  in  Figure  14  provided  the  necessary  data  to  model  the  signal 
in  PSPICE  and  determine  the  magnetic  coupling  loss  of  the  circuit.  The  peak  voltage 
output  was  16.4  millivolts  in  both  figures.  The  signal  was  a  dampened  sinusoid.  The 
following  discussion  compares  the  PSPICE  simulation  to  the  measured  calibration 
output. 
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Figure  14 


Oscilloscope  Screen  Capture  of  Input/Output  Waveforms  (Test  Coil 
Inside  Pressure  Chamber). 
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Figure  15 


PSPICE  Model  of  Calibration  Signal  (Test  Coil  Inside  Pressure 
Chamber). 


Knowing  the  resistanee  of  the  internal  eoil  was  R  =  0.02  Q ,  the  following  equation 
was  used  to  determine  the  PSPICE  parameters, 


I{t)  =  ^  , 


(3.1) 
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where: 


•  I{t)  is  the  output  eurrent  with  respeet  to  time, 

•  I ^  is  the  initial  eurrent  at  time  zero, 

•  is  the  effeetive  resistanee  of  the  cireuit, 

•  L  is  the  induetanee  of  the  test  eoil,  and 

•  t  is  the  time. 


The  first  step  determined  the  eapaeitanee,  C .  This  was  aeeomplished  by  using  the 
resonant  frequeney,  /  =  42.3  kHz  ,  from  Figure  14  and  applying  to  the  resonant 
frequeney  equation  to  determine  the  eapaeitanee. 


/  = 


>C  = 


iTTyflc  {27rfy  L 


=  0.26pF . 


(3.2) 


The  seeond  step  determined  the  eurrent,  ,  at  the  first  voltage  peak  in  Figure  14, 
where  =16.4  millivolts. 


/ 


pk-l 


V 

^^  =  0.82  A. 
R 


(3.3) 


The  third  step  determined  the  eurrent  at  the  next  peak,  , 


/ 


pk-2 


V ,  , 

^^  =  0.55  A. 
R 


(3.4) 


The  time  differenee  between  peaks  was  =  1 1 .22  ps .  Equation  (3.1)  ean  be 

rearranged  to  determine  the  effeetive  resistanee  of  the  eombined  two-eoil  pair, 

R^  =  2491.7  .  The  resistor,  RA  ,  shown  in  Figure  12  was  determined  by  adjusting  its 

value  to  make  the  output  of  the  primary  eoil  equal  to  the  experimental  value  of  8.5  volts. 

The  only  variable  unknown  in  the  PSPICE  model  was  the  eoupling  faetor.  The 
eoupling  faetor  was  iteratively  solved  for  using  the  PSPICE  model.  This  was  done  by 
eomparing  expeeted  output  of  Eigure  14  to  the  PSPICE  model  and  resetting  the  eoupling 
faetor  until  a  mateh  of  signals  was  aehieved.  The  coupling  factor,  ,  was  found  to  be 

0.00018  .  The  importance  of  the  coupling  factor  is  based  upon  the  test  coil  configuration. 
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The  test  eoil  had  to  be  inserted  through  a  one-quarter  ineh  diameter  opening  on  the  end  of 
the  pressure  ehamber.  The  geometry  of  the  test  eoil  was  a  reetangle  with  a  length  of 
seven  inehes  and  width  of  y4  of  an  ineh.  (The  resulting  area,  ,  was  5.25  square 
inehes.)  The  sample  ehamber  had  an  area,  ,  of  13.6  square  inehes.  The  fuel 

eelTs  magnetie  eross-seetion  area,  A^ ,  was  16  square  inehes.  The  expeeted  eoupling 
faetor  for  the  test  eoil,  ,  was  found  by  dividing  the  smaller  area  by  the  larger  area. 


=  ■ 


A^ 


=  0.328. 


(3.5) 


The  expeeted  eoupling  faetor  between  the  sample  ehamber  and  fuel  eell  was,  =  0.85 . 

The  large  differenee  between  the  theoretieal  eoupling  faetor  and  that  found  in 
PSPICE  led  to  a  seeond  ealibration  test  with  the  eoil  wound  two  times  around  the  aetual 
fuel  eell  windings.  This  should  have  made  the  eoupling  faetor  equal  to  one.  Figure  16 
shows  the  voltage  input  signal  at  8.5  volts  and  voltage  output  signal’s  first  peak  at  1  volt. 
The  aetual  signal  output  is  shown  in  Figure  17.  The  output  signal’s  first  peak  voltage 
was  1  volt  and  the  eoupling  faetor  from  the  PSPICE  model  was  0.012  .  This  eoupling 
faetor  was  then  used  to  normalize  the  eoupling  faetor  found  for  the  eoil  inside  the 
ehamber.  The  result  was  a  eoupling  faetor  of  0.15  .  These  results  provide  a  margin  of 
error  for  the  predieted  values  of  the  fuel  eell  experiments  with  a  range  of 
0.15  <  <  0.85  to  aeeount  for  eoupling  faetor. 
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Figure  16  PSPICE  Model  of  Calibration  Signal  (Test  Coil  Outside  Pressure 

Chamber). 
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Figure  17  Oscilloscope  Screen  Capture  of  Input/Output  Waveforms  (Test  Coil 

Outside  Pressure  Chamber). 


The  complete  lab  setup  is  shown  in  Figure  18.  The  HP  oscilloscope  is  in  the 
foreground  with  the  yellow  air  bottle  on  the  right.  The  round  metallic  cylinder  was  used 
to  diminish  the  discharge  blast,  and  the  sandbags  were  used  to  help  reduce  mechanical 
vibration.  The  signal  generator  on  the  left  of  Figure  18. 
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Figure  18  Experiment. 


3,  Operation 

A  stable  0.1  Tesla  magnetie  field  was  generated  using  four  permanent  magnets 
situated  on  eaeh  side  of  a  stainless  steel  eylinder.  A  Hall-effect  Gauss-meter  verified  the 
static  magnetic  field  in  the  gap  at  atmospheric  pressure  was  0.1  Tesla  agreeing  with  the 
calculation  and  modeling  (the  modeling  predicted  0.1  Tesla).  A  22-gauge  wire  coil  with 
515  turns  was  placed  around  the  sample  cylinder  to  sense  magnetic  flux  changes.  Air  at 
1800  psi  was  pulsed  into  the  sample  cylinder  using  valve  two.  The  gas  enhances  the  lines 
of  flux  within  the  gap  causing  a  voltage  to  be  induced  in  the  coil.  The  gas  then  exits  the 
cylinder  causing  the  lines  of  flux  in  the  gap  to  decrease.  This  flux  decrease  causes  an 
opposite  polarity  voltage  to  be  induced  in  the  coils.  The  coil  is  connected  to  a  HP 
INFINIUM  400  MHz  oscilloscope  (Figure  19),  which  measured  the  output  voltage. 
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Figure  19  HP  INFINIUM  400  MHz  Oscilloscope. 

The  microphone  connected  to  the  oscilloscope  trigger  function,  activated  upon 
receiving  an  un-amplified  signal  of  10  millivolts.  The  noise  burst  lasts  for  approximately 
200  milliseconds. 

C.  EXPERIMENTAL  RESULTS 

The  output  of  the  magnetic  fuel  cell  is  shown  in  Figure  20.  The  peak  voltage  was 
6.3  millivolts  and  lasted  for  0.3  milliseconds.  The  expected  voltage  was  2.8  millivolts  for 
a  duration  of  3.6  milliseconds  (from  Appendix  D).  The  measured  voltage  was  2.3  times 
greater  than  expected  while  the  signal  duration  was  one  tenth  of  the  expected  value. 
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Figure  20  HP  Oscilloscope  Screen  Capture  of  Magnetic  Fuel  Cell  Output  using 

Compressed  Air. 


Figure  21  shows  the  theoretical  voltage  values  versus  the  predicted  gas  output  time  (blue 
line).  The  predicted  voltage  spans  three  orders  of  magnitude  over  ten  milliseconds.  The 
theoretical  derivation  showed  a  discharge  time  of  3.6  milliseconds  (red  square).  The 
manufacturer  of  the  one-quarter-inch  valves  [11]  established  a  time  of  0.5  millisecond  for 
the  gas  to  exit  the  pressure  chamber  at  1800  psi  (blue  triangle  with  the  voltage  set  to  the 
measured  value  of  6.3  millivolts).  The  actual  voltage  and  time  were  6.3  millivolts  and 
0.3  milliseconds  (green  circle).  The  experimental  results  demonstrate  the  challenge  in  the 
analysis.  A  small  difference  in  the  expected  time  for  the  gas  to  exit  causes  the  output 
voltage  to  change  considerably. 
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Figure  21  Voltage  Output  vs.  Change  in  Time  Out. 

The  predicted  voltage  varies  linearly  with  pressure  as  shown  in  Figure  22  (blue  line). 
The  voltage  output  with  respect  to  pressure  was  very  stable  and  easy  to  predict.  The 
measured  value  (green  triangle)  was  within  a  factor  of  two  of  the  predicted  value  (red 
square). 


■  Predicted  Value 
•  Measured  Value 
Manufacturer  Time 


Figure  22  Plot  of  Predicted  Voltage  vs.  Supply  Pressure. 
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One  challenge  that  became  clear  early  on  in  the  experimentation  was  properly 
triggering  the  scope.  The  vibration  audio  signal  created  by  the  rapid  depressurization  of 
the  air  in  the  pressure  chamber  was  used  to  trigger  the  scope.  The  vibration  additionally 
induced  a  voltage  signal  from  movement  of  the  coil.  The  rapid  expansion  of  the  gas  out 
of  the  chamber  caused  the  magnets  to  vibrate  for  a  period  of  200  milliseconds.  A 
reduction  of  the  vibration  noise  magnitude  was  finally  achieved  by  tightly  wrapping  the 
magnetic  fuel  cell  with  tygon  tubing  and  placing  several  20-pound  sandbags  on  top  of  it. 


Bernoulli’s  principle  was  adapted  from  its  true  purpose  of  modeling  fluid  flows  to 
determine  timing  of  the  gas,  .  This  adaptation  and  the  mechanical  vibration 


prevented  an  accurate  measurement  of  .  The  initial  experiments  were  obtaining 

approximately  15 -millivolt  vibration  pulse  signals  shown  in  Figure  23.  This  was  later 
reduced  with  tygon  tubing  and  sand  bag  modifications  mentioned  above. 


Figure  23 


TIP  Oscilloscope  Screen  Capture  of  Vibration  Caused  Noise  Output. 


The  next  section  discusses  the  knowledge  learned  from  this  experimental  setup. 
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D.  LESSONS  LEARNED 


The  next  experimental  model  should  incorporate  the  following  improvements  to 
its  configuration. 

•  Reduce  the  magnetic  gap  and  height  to  optimize  magnetic  field  strength  in  the 
gap  and  reduce  the  air  volume  required  to  change  the  magnetic  field. 

•  Use  a  high  frequency  soft  magnetic  core  as  the  magnetic  return  path  to 
decrease  magnetic  coupling  losses  and  allow  the  coil  to  be  moved  from  the 
gap  to  the  core. 

•  Use  10-14  gauge  wires  as  the  windings  to  optimize  the  power  factor  term  and 
thus  optimize  fuel  cell  power  output. 

•  Electro-magnetically  isolate  the  entire  fuel  cell  system  and  permanently  pot 
the  fuel  cell  to  minimize  mechanical  vibration. 

•  A  software-modeling  program  such  as  FEMLAB  for  MATEAB  [13]  may  be 
useful.  This  program  allows  multi-physics  modeling  and  is  very  intuitive. 

•  Use  flexible  pressure  hoses  vice  copper  tubing  as  input  to  the  fuel  cell  to 
reduce  mechanical  vibration. 

•  Keep  the  size  of  the  valves  and  piping  larger  than  the  cross-sectional  area  of 
the  pressure  chamber  to  maximize  gas  flow. 

•  Find  solenoid-operated  high-pressure  valves  to  act  as  inlet  and  outlet  to  the 
fuel  cell.  (This  will  allow  better  signal  isolation  and  better  characterizing  of 
the  fuel  cell.) 

•  Re-examine  Bernoulli’s  Principle  and  its  application  to  the  fuel  cell  to  better 
isolate  timing  performance. 

•  Purchase  a  flow  detector  for  the  fuel  cell. 
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Chapter  IV  discussed  the  setup,  calibration  and  testing  the  magnetic  fuel  cell. 
Analysis  of  the  test  results  showed  a  close  agreement  between  the  predicted  and 
measured  voltages  and  gas  flow  times.  The  final  section  of  this  chapter  provided  an 
extensive  list  of  lessons  learned  and  suggestions  in  designing  the  next  magnetic  fuel  cell 
experiment.  Chapter  V  provides  a  conclusion  to  the  thesis,  discussing  what  a  next 
generation  magnetic  fuel  cell  is  capable  of  and  provides  additional  ramifications  of  the 
theory  behind  the  magnetic  fuel  cell. 
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V.  CONCLUSION,  OPTIMAL  DESIGN  AND  RAMIFICATIONS 


A.  CONCLUSION 

This  thesis  discussed  the  development,  design  and  testing  of  a  magnetic  fuel  cell. 
The  experiment  showed  the  magnetic  fuel  cell  produced  a  significant  voltage  output.  The 
experiment  provided  insight  to  design  a  more  advanced  experiment.  Two  lessons  learned 
were  the  importance  of  matching  the  magnetic  circuit  and  the  electrical  circuit.  An  im¬ 
properly  matched  magnetic  circuit  produced  coupling  losses.  An  improperly  matched 
electrical  circuit  did  not  maximize  the  power  transfer  of  the  magnetic  fuel  cell.  A  future 
experiment  should  be  designed  to  develop  a  larger  power  output  and  reduce  the  coupling 
and  mechanical  problems  found  in  the  first  experiment.  The  next  section  will  discuss  the 
next  generation  design  of  the  magnetic  fuel  cell. 


B,  PROPOSED  NEXT  GENERATION  DESIGN 

Figure  24  shows  the  voltage  output  of  the  ideal  magnetic  fuel  cell  using  solenoid 
operated  charge  and  discharge  valves  operating  at  high  frequency.  The  theoretical 
equations  developed  in  Chapters  If  and  Iff  were  used  to  design  an  optimized  model.  A 
10,000-psi  air  supply  tank  with  a  storage  capacity  of  220  liters  acted  as  the  gas  supply. 
The  geometry  of  the  magnetic  fuel  cell  was  a  height  of  one  millimeter  and  width  of  50 
millimeters  using  a  soft  magnetic  core  return  path.  The  gap  was  0.5  millimeters  using 
two  Nd-B-Fe-40  permanent  magnets.  The  soft  magnetic  core  used  a  20-gauge  coil  with 
10,000  windings.  The  voltage  signal  with  a  peak  of  ±  0.1 19  volts  is  a  sinusoid  operating 
at  1.6  Hz. 
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Figure  24  Voltage  vs.  Time  for  the  Optimal  Fuel  Cell. 

Figure  25  shows  the  fuel  eell  power  versus  duration  time  of  the  supply  tank.  This  plot 
superimposes  three  values  of  gap  width  varying  from  0.1  mm  to  10  mm.  To  maximize 
power  output  the  gap  width  should  be  made  as  large  as  possible. 


60 

^  50 

(A 

re  40 

^30 

o 

^  20 

O 

Q-  10 

0 


- -t— 

1 

1 

1 

♦  Gap  Width-0.1  mm 
—■—Gap  Widths  1  mm 

A  Gap  Width=10  mm 

1 

1 

1 

% 

% 

% 

% 

\ 

\ 

X 

X 

s 

. 

. . . 

— - 

10 


100 


1000 


10000 


Time  (days) 


Figure  25  Optimal  Fuel  Cell;  Power  vs.  Total  Time  (10,000  psi  220L  tank). 


Figure  26  shows  effieieney  as  a  funetion  of  fuel  eell  gap  width  for  the  fuel  eell.  Three 
gap  widths  are  superimposed  on  the  plot.  A  gap  width  of  0. 1  mm  produees  the  largest 
effieieney.  This  is  in  eontrast  to  Figure  25  whieh  showed  a  larger  gap  width  is  better.  To 
optimize  for  both  effieieney  and  power,  a  eonliguration  having  several  fuel  eells  staeked 
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together  would  be  optimal.  This  technique  is  already  employed  in  hydrogen  fuel  cell 
systems  [17]. 
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Figure  26  Optimal  Fuel  Cell:  Efficiency  vs.  Power  (10,000  psi  220L  tank). 


C.  RAMIFICATIONS 

The  theory  developed  in  this  thesis  indicates  several  areas  of  research.  Using  the 

equations  derived  in  Chapters  II  and  III,  it  has  become  clear  that  it  is  feasible  to  construct 

an  electromagnetic  gas-flow  control  valve.  Reducing  the  length  of  the  gap  of  a  magnetic 

fuel  cell  will  increase  the  magnetic  differential  pressure  (Equation  (2.44)).  With  this 

configuration,  it  should  be  possible  to  raise  the  magnetic  differential  pressure  higher  than 

the  gas  pressure  electronically;  stopping  gas  flow.  By  reducing  the  magnetic  differential 

pressure  less  than  that  of  the  gas  pressure,  the  gas  would  flow  at  a  reduced  rate.  This 

effect  would  throttle  the  gas  flow.  A  throttle  valve  using  a  mechanical  interface  to  slow 

or  stop  gas  flow  performs  this  function.  Proper  study  of  the  magnetic  pressure  effect  may 

provide  a  device  that  does  not  require  a  mechanical  interface  to  control  gas  flow.  A  non- 

contact  electromagnetic  valve  would  be  ideal  in  many  applications  where  a  mechanical 

interface  is  undesired  or  not  feasible  (e.g.,  high-pressure  high-speed  impulse  jets).  If  this 

technology  could  be  developed,  it  would  have  many  applications  both  in  the  military  and 

commercial  industry.  Such  a  electromagnetic  valve  would  benefit  the  magnetic  fuel  cell 

development  by  replacing  the  solenoid-operated  input  and  output  valves.  A  similar 

approach  has  been  realized  in  the  form  of  a  magnetic  wind  blower  built  by  a  group  of 

Japanese  researchers.  They  have  shown  that  it  is  possible  to  produce  wind  using  the 
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changing  susceptibility  of  air  and  a  high  strength  electromagnets  in  their  research  lab 

[14]. 


One  avenue  of  optimization  not  discussed  was  that  of  optimizing  the  area  term  of 
Faraday’s  Induction  Equation  (Equation  (1.24)).  Since  a  sphere  is  the  largest  volume  for 
a  fixed  surface  area  for  any  three-dimensional  object,  the  magnetic  fuel  cell’s  optimal 
shape  should  be  spherical.  Area  can  be  maximized  while  minimizing  the  number  of 
windings,  change  in  magnetic  field  and  change  in  time. 

The  Earth  is  an  example  of  a  very  large  magnetic  fuel  cell.  The  voltage  difference 
created  by  this  magnetic  fuel  cell  would  be  between  the  Earth’s  ground  and  the 
ionosphere.  Richard  Eeynman  discussed  the  mystery  behind  the  energy  in  the  ionosphere 
in  his  Eeynman  Eectures  on  Physics  [15], 

Although  the  electric  current-density  in  the  air  is  only  a  few 
micromicroamperes  per  square  meter,  there  are  very  many  square  meters 
in  the  earth’s  surfaee.  The  total  electric  current  reaeh  the  earth’s  surface  at 
any  time  is  very  nearly  constant  at  1800  amperes.  This  current,  of  course, 
is  “positive”  -it  carries  plus  charges  to  earth.  So  we  have  a  voltage  supply 
of  400,000  volts  with  a  eurrent  of  1800  amperes  -  a  power  of  700 
megawatts! 

With  such  a  large  current  coming  down,  the  negative  charge  on  earth 
should  soon  be  discharged.  In  fact,  it  should  take  only  a  half  an  hour  to 
discharge  the  entire  earth.  But  the  atmospheric  electric  field  has  already 
lasted  more  than  a  half-hour  since  its  discovery.  How  is  it  maintained? 

What  maintains  the  voltage?  And  between  what  and  Earth?  There  are 
many  questions.  [15] 

The  earth’s  magnetic  field  acts  as  the  magnetic  source  in  this  magnetic  fuel  cell 
model.  The  changes  in  air  pressure  caused  by  either  heating  or  cooling  in  the  earth’s 
atmosphere  would  cause  a  voltage  between  the  ionosphere  and  earth  (Equation  (1.37)). 
Alternatively,  depending  if  the  change  in  pressure  is  vertical  or  lateral,  a  potential 
difference  between  two  laterally  different  air  masses  would  be  produced.  Thus,  the 
implication  of  the  earth  as  a  fuel  cell  provides  a  possible  explanation  of  the  origin  of  the 
electrical  energy  in  the  earth’s  ionosphere. 

Einally,  the  magnetic  fuel  cell  when  optimized  should  provide  a  safe,  inexpensive 
alternative  to  the  hydrogen  fuel  cell  and  could  be  a  replacement  to  batteries.  Application 
of  the  magnetic  fuel  cell  concept  to  energy  storage  and  non-invasive  gas  flow  detection 
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are  two  of  many  possible  future  developments.  The  magnetie  fuel  eell  provides  a  safe, 
inexpensive,  environmentally  sound  alternative  to  traditional  portable  energy  systems  and 
will  aet  as  a  eomparison  benchmark  for  existent  and  future  energy  systems. 
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APPENDIX  A:  FUEL  CELLS 


A  fuel  cell  is  an  advanced  power  generation  device  developed  to  reduce  the 
logistics,  weight,  cost,  and  size  of  traditional  power  sources.  Alexander  Grove  invented 
an  electric  fuel  cell  in  1839  [16].  An  electric  fuel  cell  is  a  device  that  converts  chemical 
energy  into  electrical  energy  directly.  A  hydrogen  gas  is  directed  down  a  proton  ex¬ 
change  membrane,  at  the  same  time  air  flows  in  parallel  to  hydrogen  down  an  opposing 
membrane.  In  between  the  two  membranes,  an  electrolyte  transports  the  proton  given  up 
by  the  hydrogen  to  the  oxygen.  Energy  is  produced.  The  waste  product  of  this  process  is 
water.  Similar  to  a  battery,  each  fuel  cell  stack  is  electrically  connected  in  series  to 
achieve  a  higher  voltage.  Unlike  chemical-reaction  batteries,  fuel  cells  last  as  long  as  the 
supply  of  oxygen  and  hydrogen  lasts. 

Electric  fuel  cells  are  used  for  a  variety  of  power  applications  to  include:  cell 
phones,  cars,  bikes,  submarines,  UAVs,  satellites,  commercial  and  residential  power 
production.  A  Proton  Exchange  Membrane  (PEM)  type  Ballard  Euel  Cells  fuel  cell  was 
employed  successfully  in  the  late  1980’s  by  Canada  in  a  commercial  submarine  [17]. 
Congress,  within  the  National  Defense  Authorization  Act  for  fiscal  year  1998,  requested 
DoD  to  report  on  the  development  and  application  of  fuel  cells  [2].  DoD  has  established 
the  DoD  Euel  Cell  Demonstration  Program.  This  program  is  intended  to  stimulate 
growth  and  economies  in  the  fuel  cell  industry.  DoD  also  established  the  Euel  Cell  Test 
and  Evaluation  Center  run  by  the  US  Army  Corps  of  Engineers.  This  center  is 
responsible  for  validating  fuel  cell  technology. 

There  have  been  several  initiatives  in  the  last  decade  to  determine  the  viability  of 
fuel  cell  usage  within  the  DoD.  These  include:  US  Coast  Guard  test  of  Molten  Carbonate 
fuel  cell  technology  aboard  the  CGC  Vindicator  in  1994  [18];  ONR-sponsored  ship 
impact  studies;  a  500  Watt  sulfur  solid  oxide  fuel  cell  technology  test  in  1997  started 
through  a  Navy  SBIR;  and  a  Marine  Corps  deployment  of  light  weight  Proton  Exchange 
Membrane  Technology  in  exercises  in  Hawaii  in  1999  and  2000  [19].  There  are  over  50 
DoD  installations  now  using  fuel  cells  [2].  Most  initiatives  have  shown  favorable  results 
using  the  fuel  cells  including  indication  of  lower  fueling  costs  for  fixed  land  sites  and 
increased  ranges  for  vehicles  using  fuel  cells. 
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APPENDIX  B  DERIVATION  OF  COUNTER  EMF 


This  Appendix  shows  the  basie  derivation  of  the  eounter  EMF  produeed  by  the 
magnetie  fuel  eell.  The  eross-seetion  of  eore  is  shown  in  Figure  27.  The  equation  for 
magnetie  flux  density  produeed  from  a  straight  wire  of  length,  L  ,  and  direet  eurrent,  /  , 
is  shown  by  [10] 


B  = 


L 


+  r 


(4.1) 


where  r  is  the  perpendieular  distanee  from  the  loop  edge  to  the  eenter  of  the  loop. 
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Figure  27  Cross  Seetion  of  the  Soft  Magnetie  Core. 
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The  magnetie  flux  density  from  a  reetangular  loop  is  the  sum  of  the  magnetie  flux 
densities  of  the  four  sides  developed  using  Equation  (4.1).  Therefore,  for  side  ,  the 
magnetie  flux  density  beeomes 
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The  magnetic  flux  density  for  side  L  becomes; 


B  = 


^n- 


f^e^ 

2 

c 

[2j 

[2) 

(4.3) 


Due  to  symmetry  the  magnetic  flux  density  for  each  parallel  wire  can  be  doubled 
resulting  in  the  flux  density  at  the  center  of  the  center  ’  Equation  (4.4)  or  re¬ 
arranged  for  current,  Equation  (4.5), 
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APPENDIX  C:  DERIVATION  OF  THE  PREDICTED  VALUES 


This  Appendix  determines  the  predicted  output  parameters  of  the  magnetic  fuel 
cell  experiment  using  equations  derived  in  Chapter  III.  Three  steps  are  performed.  Step 
one  determines  the  timing  of  the  magnetic  fuel  cell.  Step  two  determines  the  magnetic 
parameters  of  the  magnetic  fuel  cell  using  magnetic  circuit  analysis.  Step  three  uses  the 
parameters  found  in  the  previous  two  steps  to  determine  the  predicted  output. 

The  first  step  is  to  determine  the  initial  parameters  in  the  gas  chamber  (note; 
actual  values  and  unit  have  been  provided  to  aid  understanding  of  the  development  and 
ease  unit  analysis): 


Ps  =  P,a,  =  34.5  MPa 
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= - M  = 
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^3391  =  1.198^ 
m 


m  =  n  V  = 

atm  r  atm  gap 


m  —  n  Y  — 

gap  r'  gap  gap 


kg 


1.198^ 
mV 


1.5 X 10'"  m^=  1.8x10  "kg 


kg 


410.7^ 

V  mV 


1.5x10'"  m'  =0.0616  kg. 


The  next  step  is  to  determine  the  time  for  the  gas  to  exit  the  chamber. 


1 — 

V  — 

^  P 

g^P 

P  ^ 

atm 

9 

34.5  MPa 

101325  Pa 

-35- 

^  valve!  -il 

{Pgapl 

P atm  J 

410.7  3^ 

1.198  3^ 

s 

V  m' 

m'  ^ 

=  410.7  i|(0.001 171  -  =  16.8^2 

m  ^  3  s  s 


m 

At 


m„ 


0.0616  kg 

16.8  tS 


=  0.00366  s 
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The  timing  of  the  fuel  cell  has  been  determined;  the  next  series  of  steps  analyzes 
the  magnetic  circuit  portion  of  the  magnetic  fuel  cell.  This  analysis  starts  off  with  the 
determination  of  the  operating  point  of  the  magnetic  circuit  [5]. 

First  the  susceptibility,  %  >  is  determined  and  the  permeability,  // ,  of  the 
experiment. 


^  0.0508 

~  =  3.14  X - =  0.08  m 


/?,  =/?„,=  0.2032  m 
w  =  w  =  0.0508  m 


X  gap 


410.7x0.2 


M 

m'; 


3449x10^ 


gap 


0.032 


kg 

mol 


-^10^-^^  =  0.000111 


cm 


=  i  +  =  1-00011 


( 


Xo. 


kg 


1.198^x0.2 


m 


3449x10 


-6  cm 


0.032 


kg 

mol 


-^^10 -^^  =  3.25x10  ’ 


cm 


/^..  =  1  + Jo  =1-0000003245. 


Next  determine  the  magnetic  flux  density  into  the  gap,  B.^ , 


^  11600  A  , 

= - =  922832  —  conversion ) 

"  0.01257  m  m  ^  ’ 

nwX  (0.0508x0.2032)  _  ^  _ 

2Peoreslh[^  ^  Vo.2032 V  0.0508^ 


H 


m 


-B 


r 


PcLm 

^gap  _|_  ^core 
Pgap  P  air 


^loop 


+ 


A 

Ha 
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H.= 


-1.25  T 


4;z-xl0~'  X  0.0127  m  1.25  T 

0.0508  m  ^  0.08  m  ^  922832  — 

1.00011  1.000000325  m 


//  =-846571  — 

m 

m 


A 


4;rxl0  x  846571  — x  0.0127  m 


= 


m 


0.0508  m  ^  0.16  m  0  077 

1.00011  ^1.000000361” 


S=  0.100981  T. 


Second,  determine  the  magnetic  flux  density  as  the  gas  chamber  is  discharged. 
A,.,,  =l  +  ;ro,  =1-000000325 
H^  = 


-B.. 


BoK 


L 


L 


Ba 


-1.25  T 


4;r  X 10  x  0.0127  m 


0.0508  m 


0.08  m 


+ 


1.25  T 


A 

922832  — 
1.000000345  1.000000345  m 


//  =  -846574  — 

m 

m 


A 


4;rxl0  x  846574  —  x  0.0127  m 


D  _  _ 

0.0508  m 


m 


0.08  m 


1.000000345  1.000000345 

5  =0.100977  T 

A5  =  5„„,-5,„=4hT. 


-.077  m 
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Therefore  the  measured  voltage  ean  be  ealeulated,  along  with  the  eurrent  and 
induetanee  [20]; 


V,=-n  —  A  =  -515x  X  0.0102  m"  =-5.5  mV 


R  = 


At,.„  0.00366  s 

=  31.99 


1  1 

+ 


32  Q  10^  Q 

r.  -5.5  mV  X  31.99 

—ed  32f2  +  31.99f2 

J  ^  Vgut-V—ed  ^  -5.5  mV +  2.75  mV  ^  ^ 

R  32  Q. 


L  =  ^^a 


f8a^ 

In 

-1.75 

Vbj 

T  =  4;r  X 10  x  0.05 


In 


8x0.05 

0.0002 


-1.75 


:(515)"  =97.5mH 
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APPENDIX  D:  DERIVED  PARAMETERS 


Appendix  D  shows  the  eomplete  list  of  input  data,  eonstants  and  derived 
parameters  for  the  experiment. 


Table  2.  Table  of  Derived  Input  Values  and  Constants. 


Input  Values 

Magnitude 

Units 

Supply  Gas  Pressure 

12,400,000 

Pa 

Supply  Gas  Temperature 

293 

K 

Atmospheric  Pressure 

101,352 

Pa 

Gap  Volume 

150 

m= 

Width  of  Magnets 

0.0508 

m 

Height  of  Magnets 

0.2032 

m 

Length  of  Magnets  (combined) 

0.0125 

m 

Length  of  Gap 

0.0508 

m 

Coercivity 

922,832 

A/m 

Residual  Magnetism 

1.25 

T 

Number  of  Windings 

515 

- 

Radius  of  Wire 

0.0002 

m 

Constants 

Ideal  gas  constant 

8.314 

J/(mol-K) 

Molar  Mass  of  Nitrogen  Gas 

0.028 

kg/ mol 

Molar  Mass  of  Oxygen  Gas 

0.032 

kg/mol 

Experimental  Susceptibility  of  Oxygen 

0.0034 

Permeability  of  Free  Space 

0.00000126 

H/m 
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Table  3.  Table  of  Derived  Output  Values. 


Output  Parameters 

Predicted 

Units 

Supply  Tank  Volume  (at  1800  psi) 

0.37 

m3 

Molar  Mass  of  Air 

0.029 

kg/mol 

Mass  of  Gas  in  Supply  Tank 

152 

kg 

Density  of  Gas  in  Supply  Tank 

410.7 

Density  of  atmosphere 

1.198 

Mass  of  Gas  in  Gap  (at  14.7  psi) 

18,000 

kg 

Mass  of  Gas  in  Gap  (at  1800  psi) 

0.0616 

kg 

Density  of  Gas  in  Gap  (at  1800  psi) 

410.7 

kg/m^ 

Pressure  of  Gap  in  Gap 

12,400,000 

Pa 

Velocity  of  the  Gas  Entering  Valve  2 

35 

m/s 

Area  Gap  Gas 

0.001171 

m^ 

Mass  Flow  Rate  of  Gas  Out  of  Gap 

16.8 

kg/s 

Time  for  Gas  to  Discharge 

0.00366 

s 

Length  of  Magnetic  Return  Path 

0.08 

m 

Height  of  Magnetic  Return  Path 

0.2032 

m 

Width  of  Magnetic  Return  Path 

0.0508 

m 

Density  of  Oxygen  in  Gap  (at  1800  psi) 

82.14 

Susceptibility  of  Gap  (at  1800  psi) 

0.000111 

- 

Permeability  of  Gap  (at  1800  psi) 

1.00011 

- 

Density  of  Oxygen  (at  14.7  psi) 

0.2396 

Susceptibility  of  Gap  (at  14.7  psi) 

0.00000033 

- 

Permeability  of  Gap  (at  14.7  psi) 

1.00000032 

- 

Kloop 

0.077 

m 

Magnetic  Field  Strength  of  Magnets  (1800  psi) 

-846571 

A/m 

Magnetic  Flux  Density  (1800  psi) 

-0.100981 

T 

Magnetic  Field  Strength  of  Magnets  (14.7  psi) 

-846574 

A/m 

Magnetic  Flux  Density  (14.7  psi) 

-0.100977 

T 

Change  in  Magnetic  Flux  Density 

4 

uT 

Area  Gap  Magnetic 

0.0102 

m^ 

Voltage  Output  of  Magnetic  Fuel  Cell 

-5.5 

mV 

Effective  Coil  Radius 

0.05 

m 

Effective  Load  Resistance 

32.0 

ohms 

Measured  Voltage 

-2.75 

mV 

Current  Output  of  Magnetic  Fuel  Cell 

-86 

mA 

Inductance  of  Coil 

97.5 

mH 

APPENDIX  E:  LIST  OF  SYMBOLS 


coercivity 

X 

susceptibility 

Mr 

relative  permeability 

B 

magnetie  flux  density 

H 

magnetic  field  strength 

Mo 

permeability  of  free  spaee 

K 

experimental  suseeptibility 

M 

molar  mass 

P 

material  density 

M 

magnetization 

M 

magnetie  permeability 

K 

number  of  atoms  per  unit  volume 

q 

electron  charge 

2 

<r  > 

expeetation  value  of  orbital  area 

z 

atomie  number 

eleetron  mass 

N, 

nitrogen 

Br^ 

bromine 

Bi 

bismuth 

0^ 

oxygen 

Mm 

atomie  magnetie  moment 
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Ks 

T 

3 

dL 


n 


I 

4 


L 


A 

(!> 

V 

R 

AE 

AE 

AE 


KE 


PE 


V 


/ 


g 

hf 

m 


Boltzmann  constant 

temperature 

magneto-motive  foree 

unit  length 

number  of  windings 

induced  current 

eore  magnetic  field  strength 

mean  length  of  core 

gap  magnetie  field  strength 

gap  length 

eross-sectional  area  of  magnetie  material 

magnetic  flux 

reluctance 

voltage 

resistance 

change  in  kinetic  energy  of  gas 
change  in  potential  energy  of  gas 
change  in  pressure  energy  of  gas 
final  velocity  of  gas 
initial  velocity  of  gas 
gravity 

final  height  of  gas 
mass  of  the  gas 
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Pi 

V 

AP 

R 

r 

A(j) 

At 

AB 

Xf 


Xi 


f 


H 

Hi 

h 


valvel 


Ps 


initial  height  of  gas 
gas  volume 
final  pressure  of  gas 
initial  pressure  of  gas 
gas  veloeity 

change  in  pressure  of  gas 
number  of  moles 
gas  constant 
ideal  gas  constant 
flux  change 
change  in  time 

change  in  magnetic  flux  density 
final  susceptibility 
initial  susceptibility 
change  in  susceptibility 
final  magnetic  field  strength 
initial  magnetic  field  strength 
height  of  gap 

velocity  of  gas  entering  valve  one 
velocity  of  the  gas  in  the  supply  tank 
pressure  of  the  gas  in  the  supply  tank 
density  of  gas  in  the  supply  tank 
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gap 


p. 


A. 


'valvel 


gap 


^  valveX 


gap— in 


valve! 


gap— out 


Pa 


A. 


valve! 


^  cycle 


‘'valvel 


‘'valve! 


fc 


cycle 


F 


cycle 


J 

ds 


pressure  of  the  gas  in  the  gap 
density  of  the  gas  in  the  gap 
cross-sectional  area  of  the  supply  valve 
cross-sectional  area  of  the  gap 
velocity  of  the  gas  entering  valve  one 
velocity  of  the  gas  entering  the  gap 
height  of  the  gap 

velocity  of  the  gas  exiting  valve  two 

velocity  of  the  gas  leaving  the  gap 

atmospheric  pressure 

atmospheric  density 

cross-sectional  area  of  valve  two 

velocity  of  the  gas  entering  the  atmosphere 

cycle  time  of  the  magnetic  fuel  cell 

time  to  open  and  shut  valve  one 

time  to  open  and  shut  valve  two 

cycle  frequency  of  the  magnetic  fuel  cell 

volume  flow  rate  of  the  cycle 

duration  time 

induced  current  density 

unit  area  normal  to  the  magnetic  flux 
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L 

m 


gap 

core 

^core 

^  core 

K 

Pgap 

f^core 

k 

K 

Br 

Pm 

He 

B,„ 


AE 


cycle 


R 


magnet  magnetic  field  strength 

length  of  magnets 

gap  magnetic  field  strength 

length  of  the  gap 

core  magnetic  field  strength 

length  of  the  core 

gap  magnetic  flux  density 

core  magnetic  flux  density 

magnet  magnetic  flux  density 

gap  permeability 

core  permeability 

simplification  constant 

core  height 

core  width 

residual  magnetization 
permeability  of  magnets 
magnet  coercivity 

magnetic  flux  density  at  high  pressure 
magnetic  flux  density  at  atm.  pressure 
electrical  energy  of  the  cycle 
load  resistance 
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cell 


AP 


mag 


AE 


total 


AE 


cell 


n 

R 

X 

y 

a 

b 


^  layer 


E 


'layer  1 


E 


'layer  2 


E 


'coil 


Pw 

A 

^wire 

m 

lo 

K 

E 


power  output  of  the  magnetie  fuel  eell  eaeh 
cycle 

magnetic  pressure  drop  caused  by  magnetic 
fuel  cell 

efficiency  of  the  magnetic  fuel  cell 

total  energy  of  the  gas 

total  energy  of  the  magnetic  fuel  cell 

power  factor  term 

horizontal  distance 
vertical  distance 
magnetic  core  radius 
wire  radius 

number  of  windings  on  a  layer 
total  number  of  layers 
length  of  first  layer 
length  of  second  layer 
total  length  of  coil 
resistivity  of  wire 
cross-sectional  area  of  wire 
output  current  with  respect  to  time 
initial  current  at  time  zero 
effective  resistance  of  the  circuit 
inductance  of  the  test  coil 
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t 

c 

f 

I 


pk-\ 


V, 


pk-l 


>/t-2 


V. 


pk-2 


At 


peaks 


A, 


A 

( 

K 

K., 


chamber 


r 


B 


center 


time 

capacitance  of  the  magnetic  fuel  cell  coil 

resonant  frequency  of  calibration  circuit 

output  current  at  peak  one 

output  voltage  at  peak  one 

output  current  of  next  peak 

output  voltage  of  next  peak 

time  difference  between  peaks 

area  of  test  coil 

area  of  pressure  chamber 

fuel  cell  magnetic  cross-section  area 

coupling  factor 

voltage  output  of  fuel  cell 

perpendicular  distance  from  loop  edge  to 
center  of  loop 

flux  density  at  the  center  of  gap 
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